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Preface 

What  has  been  learned  in  the  past  few  months  will,  undoubtedly, 
be  reflected  in  the  following  pages;  moreover,  so  will  changing  altitudes 
and  points  of  view.  Thus  what  was  solely  an  absorption  problem  at  the 
beginning  has  grown  somewhat  in  the  intervening  time  reflecting,  perhaps, 
increasing  awareness  of  the  problems  with  which  ionospheric  physics  is 
confronted,  which,  while  being  immediate  in  character,  are  always  difficult 
to  handle  with  the  utmost  preciseness  t  so  that  I  have  had  to  draw  from 
other  disciplines  and,  though  I  was  not  assured  of  more  precisenss,  these 
disciplines  did  help  in  a  greater  understanding  of  the  problem  and  its 
related  facts.  In  the  paper,  this  is  translated  into  certain  basis  features. 
First,  a  fairly  lenghy  examination  of  certain  aspects  of  magneto-ionic 
theory  has  been  included  which,  in  my  opinion,  is  essential  to  my  understanding 
of  the  calculation  of  the  absorption  integrals.  This  examination  forms  the 
core  of  Chapters  II  and  III.  While  I  realize  that  it  may  be  hard  going  for 
anybody  reading  this  paper,  its  inclusion  also  reflects  the  amount  of  time 
that  was  spent  on  it.  In  Chapter  IV,  particularization  is  made  to  the  more 
relevant  expressions  to  the  calculations  themselves.  The  rest  of  the 
Chapters  concern  themselves  with  results,  discussion  and  conclusions  along 
with  pertinent  recommendations.  Much  material  has  been  omitted  which  was 
related  to  the  topic  although  not  in  a  very  direct  way.  This  is  reflected 
in  the  Bibliography;  which  contains  all  the  relevant  literature  which  I 
had  occasion  to  refer  to  during  the  time  of  work.  It  is  hoped  that 
anybody  wanting  to  pursue  the  subject  further  may  use  the  Bibliography 
as  a  guide  of  where  to  look  first. 
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During  the  time  that  this  vcrk  was  being  done,  several  people 
cooperated  in  one  way  or  another  and  to  them  goes  my  appreciation!:  to 
Mr,  Richard  S.  Allen,  who,  with  his  tinfailing  enthusiasm  and  ever  present 
optimism,  made  an  optimist  out  of  me,  at  least  some  of  the  time;  to  Mrs, 
Isabel  Hussey,  who  helped  in  the  computer  input  and  output  formats  and 
some  of  the  graphing;  to  Miss,  Patricia. Kennan ,  who  did  the  programming, 
for  her  patience;  and  to  Mr.  Paul  Lyons,  who  did  some  of  the  graphs 
quickly  when  there  was  not  too  much  time.  To  all  others,  whose  names  I 
may  omit,  but  were  in  some  way  instrumental  to  the  final  work,  I  also 
express  my  gratitude  and,  very  especially,  to  all  those  who,  despite  my 
better  judgement,  convinced  me  not  to  write  this  paper  in  Spanish, 

Jiu'.o  Edgardo  Barety 
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Abstract 

The  ra ague to -ionic  theory  is  examined  in  its  relation  to  absorption 
calculations.  For  representative  profiles  the  complex  index  of  refraction 
was  computed  using  the  Appleton-Hartree  and  Sen -Wy Her  magneto-ionic 
formulas.  The  height  profiles  of  the  real  and  imaginary  part  of  the  index 
of  refraction  are  plotted  and  the  absorption  corresponding  to  day-time  and 
night-time  conditions  is  computed  and  its  frequency  dependence  determined. 

It  is  determined  that  absorption  for  general  conditions  also  follows 
closely  an  inverse  square  law.  although  there  are  deviations,  particularly 
for  the  extraordinary  ray.  It  is  determined  that  general  propagation 
conditions  do  not  affect  radically  the  possible  use  of  differential 
absorption  for  total  absorption  measurements.  A  possible  method  for  doing 
this  is  given.  The  vari&i  .on  of  the  absorption  coefficient  as  a  function 
of  parameter  N  is  also  determined  and  the  linearity  condition  &*£()$*>)  A/ 
is  shown  to  hold  under  certain  conditions  for  general  propagation  conditions* 
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I.  INTRODUCTION 

In  radio  communications  and  atmospheric  sciences  in  general* 
there  has  always  bean  interest  in  determining  the  absorption  of 
radio  waves  on  traversal  through  the  ionosphere.  As  early  as  193&, 
Appleton  (Ref  4)  made  calculations  under  simplified  conditions. 

These  simplified  conditions  included  the  quasi-longitudinal 
approximation,  which  is  equivalent  to  propagation  of  the  wave  in  a 
direction  nearly  parallel  to  the  magnetic  field;  non-deviate  absorption, 
i.e.  the  index  of  refraction  M  is  nearly  1  and  a  Chapraan-I*yer  model 
for  the  absorbing  layer.  At  that  time  the  Chapman-Layer  model  was 
considered  to  be  valid  for  the  E- layer.  Later  (1946),  Jaeger  (Ref  27) 
performed  calculations  of  the  same  type  extending  Appleton's  work  to 
deviative  absorption.  He  also  considered  reflections  from  the  top 
side  for  a  Chapman-Layer.  The  absorption  for  topside  reflections 
will,  in  general,  be  different  from  that  for  reflections  from  the 
bottom-side  because  of  the  asymmetry  of  the  Chapman-Layer  around  the 
height  of  maximum  ionization.  Jaeger  retained  the  other  restrictions 
that  Appleton  used.  Another  detail  which  should  be  noted  is  that 
both  were  based  on  the  classical  magneto-ionic  formula  for  the  index 
of  refraction  derived  by  Appleton  (Ref  3)* 

The  ionosphere,  however,  is  much  complicated  than  the  use  of 
some  of  tne  previous  restrictions  would  seem  to  imply.  In  particular, 
the  electron  density  profiles  of  the  ionosphere  are  quite  different 
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from  the  Chapman-'  yyer  model.  Also  the  collision  frequency  proi'.les 
are  more  complicated  than  the  simple  exponential  decrease  used  in 
the  calculations  mentioned  in  the  first  paragraph.  The  use  of  more 
realistic  profiles  leads  to  the  use  of  electron  density  profiles 
and  collision  frequency  profiles  for  which  no  analytical  expression 
exists.  This  means  that,  in  general,  analytical  integration  of  the 
absorption  integrals  is  not  possible.  Finally,  the  quasi -longitudinal 
approximation,  though  quite  good  for  a  great  range  of  conditions 
(Ref  45 : 79)  is  only  a  particular  case  of  general  propagation  of  a~ 
wave  in  a  plasma  medium  upon  which  a  constant  magnetic  field  has 
been  imposed. 

The  author  has  endeavored  in  this  paper  to  calculate  the  absorption 
integrals  for  certain  representative  conditions  as  represented  by 
typical  electron  density  profiles  and  collision  frequency  profiles.  The 
conditions  imposed  on  the  integral  were  fairly  relaxed  in  the  general 
direction  indicated  by  the  previous  paragraph.  They  were  computed  for 
both  magneto-ionic  components  of  the  wave. 

The  General  Problem 

A  small  problem  such  as  this  one  usually  forms  part  of  a  larger 
one.  This  larger  problem  has  motivated  the  author  in  the  work  related 
to  the  present,  more  specific  problem. 

The  general  problem  is  the  determination  of  the  energy  spectrum 
of  cosmic  radio  sources  (Ref  11),  which  may  help  determine  the  energy 
mechanisms  by  which  they  emit  radio  waves.  It  is  observed  that  as 


2 


sp/ph/65-1 


the  frequency  decreases  from  very  high  values  towards  the  HF  region, 
the  spectral  intensity  increases.  It  may  be  assumed,  however, 
that  there  should  be  a  drop  in  intensity  as  the  frequency  approaches 
zero,  in  accordance  with  the  physical  laws  that  govern  the  behaviour 
of  radiating  sources.  The  shape  of  the  spectral  curves  between 
5  and  30  ivC/S  is  very  much  in  doobt  at  the  moment.  A  knowledge 
of  the  shape  of  the  energy  curves  in  the  5-3°  MS/ S  range  could  very 
well  help  determine  what  type  of  emission  mechanism  is  present; 
whether  thermal  sources  with  a  spectral  response  similar  to  that  of 
a  blackbody,  or  non-thermal  sources  like  relativistic  particles 
emitting  radiation  when  trapped  in  magnetic  fields,  are  responsible 
for  the  emission.  The  reason  why  the  behaviour  of  the  spectrum  of  the 
radio  sources  is  in  doubt  between  5  a*1**  3°  MC/S  is  that  it  is  masked 
by  the  day  to  day  variations  of  strong  absorption  in  the  lower  parts 
of  the  ionosphere .  It  is  hoped  that  an  understanding  of  the  absorp¬ 
tion  undergone  by  a  wave  passing  through  the  ionosphere  will  be  a 
positive  step  in  resolving  this  difficulty. 

The  fairly  broad  glimpse  given  here  serves  to  emphasize  the 
connection  between  the  problem  being  examined  in  this  paper  and  the 
broader  one  of  cosmic  radio  sources. 

The  Mathematical  and  Physical  Context 

Apart  from  the  manner  in  which  the  radio  spectrum  problem  fits 
into  the  overall  context  there  is  also  a  mathematical  and  physical 


3 


sp/fh/65-i 


context  on  which  the  present  calculations  are  based*  For  the  absorp¬ 
tion  calculations  presented  in  this  paper,  it  is  Decessary  to  know 
the  refractive  index  of  the  medium  and  it  can  be  calculated  from  the 
magnetc-ionic  theory.  This  theoiy  forms  the  immediate  mathematical 
and  physical  context  of  the  present  calculations.  Therefore  ,  the 
present  paper  devotes  considerable  space  to  those  aspects  of  the  theory 
which  are  essential  for  the  absorption  calculations.  The  physical 
considerations  leading  to  the  derivations  of  the  useful  equations  have 
been  examined  at  length.  These  derivations  follow  closely  those  given 
in  the  literature;  not  all  has  been  included,  however,  and  the  interested 
reader  is  referred  to  the  oroginal  sources.  Furthermore,  this  paper 
examines  the  theory  on  two  different  levels;  the  simple,  original 
point  of  view  by  Appleton  and  termed  "classical",  and  the  more 
sophisticated  approach  of  Sen  and  Wyller.  As  will  be  seen,  these 
approaches  give  the  index  of  refraction  at  a  point  as  a  function  of 
the  value  of  certain  parameters  at  that  point.  The  equations  are 
then  valid  even  for  inhomogeneous  media.  To  get  the  attenuation  of 
the  wave,  a  simplified  treatment  of  Maxwell's  equations  for 
inhomogeneous  media  exists.  It  is  called  the  WKB  method  and  the 
solutions  give  an  indication  of  the  attenuation  of  the  wave  in  a. 
slowly-varying  medium  (fief  12:  Ch.  9)»  A  few  wards  have  been  set 
down  on  the  WKB  solutions  at  the  end  of  Chapter  2. 

Therefore,  the  next  two  chapters  occupy  themselves  with  classical 
magneto-ionic  theory  and  its  generalizations.  The  fourth  chapter  is 
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about  numerical  applications  and  specifically  about  the  present 
calculations.  The  final  chapter  concerns  itself  with  a  discussion 
of  results  and  some  conclusions  and  recommendations# 
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II,  CLaSLICaL  MaGNETO-IONIC  THE  (BY 

The  magneto -ionic  theory  serves  as  the  basis  for  most  calcula¬ 
tions  in  the  ionosphere.  It  concerns  itself  with  the  propagation  of 
electromagnetic  waves  through  a  magneto-ionic  medium.  A  magneto¬ 
ionic  medium  is  defined  by  Ratcliff e  (Ref  45:^)»  as  one  in  which  free 
electrons  and  heavy  positive  ions  are  situated  in  a  uniform  magnetic 
field  and  are  distributed  with  statistical  uniformity,  so  that  there 
is  no  resultant  space  charge.  For  most  purposes  this  is  sufficient, 
although  sometimes  it  is  necessary  to  include  neutral  particles  and 
negative  ions.  To  treat  the  absorption  problem,  it  is  possible  to  use 
two  points  of  view:  The  macroscopic  and  the  microscopic.  In  the 
former  point  of  view  an  explanation  of  what  happens  to  the  elementary 
charges  comprising  the  medium  is  avoided.  Instead  tft  the  parameters 
£  ,  the  relative  capacitivity;  6  ,  the  conductivity;  and  /*»,  ,  the 
relative  permeability  are  used  to  account  for  the  effect  on  the  electric 
field  of  e^ual  numbers  of  oppositely  charged  particles  and  also  of  the 
translational  motion  of  bound  charges.  6  ,  obviously,  is  related  to 
the  translational  motion  of  the  constituent  charges.  Since  in  many 
media  these  parameters  can  be  measured  with  relative  ease,  the 
macroscopic  point  of  view  is  the  more  useful.  In  plasmas  some  modifica¬ 
tions  are,  however,  necessary  because  in  a  plasma  the  constitutive 
parameters  of  the  medium  are  not  easily  measured.  The  reason  for  this 
is  that  the  space  charges  present  and  the  boundary  condition  imposed 
by  the  measuring  apparatus  affect  each  other  significantly.  The 
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modification  then  is  to  use  a  microscopic  approach  to  determine  the 
constitutive  parameters  and  then  use  these  in  Maxwell's  equations 
to  determine  the  fields.  This  approach  is  what  mainly  constitutes 
the  Lorentz  theory. 

Classical  Magneto-Ionic  Theory 

Brandstatter  (Bef  10;40)  gives  a  brief  description  of  this  theory. 

The  Lorentz  theory  deals  directly  with  the  elementary  charges  of  the 
medium.  The  theory  dispenses  with  the  concept  of  the  medium.  Instead* 
it  visualizes  the  positive  and  negative  charges  which  constitute  the 
medium  to  be  in  free  space.  It  postulates  that  all  permeable  bodies 
are  composed  of  a  large  number  of  charged  microscopic  particles 
separated  by  free  space.  In  conducting  bodies,  these  microscopic 
particles  are  capable  of  being  moved  through  the  body  because  of  electrical 
forces;  that  is.  in  non-conducting,  or  weakly  conducting  bodies*  elastic 
forces  are  presumed  to  bind  the  particles  to  some  equilibrium  position. 

Any  displacement  that  they  may  undergo  from  their  position  of  equilibrium 
is  very  small,  and  when  displacement  occurs  the  medium  becomes  polarized* 
It  la  also  aaauuiou  that  ths  net  charge  of  the  medium  is  zero# 

Furthermore,  free  displacement  currents  are  assumed  to  exist  also 
within  the  particles  in  addition  to  their  existence  in  the  free  space 
between  the  particles.  The  way,  then,  that  the  effect  of  the  medium 
is  taken  into  account  is  to  consider  the  motions  of  the  charged 
particles  in  electromagnetic  fields.  Then  whether  the  charges  are 
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free  or  bound,  a  convection  current  can  be  definedo  If  in  a  certain 
portion  of  apace  there  are  N  particles  per  unit  volume,  each 
with  a  charge  ^  and  a  mass  >n  ,  then  if  _/  is  the  time  average 
of  the  charge  velocity,  each  will  produce  a  current  ^ jf  and  for  the 
totality  the  current  will  be  f/fY  •  Hence,  the  average  convection 
current  ^  v  ~f? : I » 

It  is  not  clear,  however,  whether  J  is  proportional  to  the 
electric  field  or  proportional  to  £  ,  and  out  of  phase  with  E  by 

(i.e.  Aether  it  is  a- conduction  current  or  a  displacement  current). 
So,  in  the  most  general  case,  it,  should  be  the  sum  of  the  free 
displacement  current  and  the  convection  current* 

Summarizing t 

The  Lorentz.  microscopic  theory  describes  the  electromagnetic 
phenomena  by  the  same  set  of  Maxwell’s  equation  as  the  macro see tic 
theory  t 


7X  E  =  -= 

(1.1) 

<7  X  H  —  J  +  ^ 

(1.2) 

v  2  r  f 

(1.3) 

7 .  B  -  0 

(1.4) 
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medium  in  this  point  of  view,  the  equations: 

(1.5) 

(1.6) 

(1.7) 

are  given  "in  vacuo", 

;  ***•  (1.8) 

■and  the  current  density  is  given  directly  in  terms  of  the  moving 
charges: 

J  =%£<%/  volume  of  region  containing  chargeo.  (1.9) 
Maxwell's  equations  relate  the  field  to  the  charges  and  their  motions, 
but  since  these  are  not  known,  it  is  necessary  to  supplement  them 

with  the  dynamical  equations  of  motion  for  the  charges.  These 

% 

equations  have  the  form,  for  each  particle,  pt  : 

where  B  is  the  total  magnetic  induction  field  evaluated  at  the 
position  ,  of*  the  particle  j  E  is  the  total  electric  field 
evaluated  at  L«  •  These  fields  can  consist  in  part  of  external 
fields,  and  in  part  of  fields  due  to  the  themselves. 

The  set  of  equations  must  be  solved  simultaneously  for  the 
dynamical  behavior  of  the  charges,  and  from  this  behavior  we  can 


but  since  there  is  no 

D  s  €  § 

J  =  6  £ 
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derive  the  constitutive  parameters  for  the  microscopic  model  of 
the  medium. 

In  the  Lorentz  theory,  the  motion  of  the  electrical  charges 

is  described  in  terms  of  the  polarization  vector  (dipole  moment  per 

unit  volume)  rather  than  in  terms  of  the  convection  current.  In 

this  case  the  polarization  vector  is  P  =  W^rr  ,  r  being  the  average 

# 

displacement  of  the  charged  particles.  Consequently,  J  =  P  . 

In  the  present  absorption  calculations  it  is  necessary  to  know 
how  the  index  of  refraction  varies  with  some  of  the  parameters  of 
the  ionosphere.  It  is  the  magneto-ionic  theory  which  provides  the 
answer. 


A.  -  THE  WAVS  EQUATION 

In  a.  magneto-ionic  medium  as  defined  by  Ratclif fe ,  the  net 
free  charge  J  =  o.  Maxwell’s  equations  then  became s 


V-B  =  0 


~  0 

sr 


vx H  -  I  + J f 


(1.3) 
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furthermore,  all  field  components  in  Maxwell's  equations  are  assumed 
to  contain  the  factor 

Uu  t - k  r ) 


A  new  notation  is  introduced  to  simplify  certain  calculations.  Vectors 

will  be  written  in  their  component  form,  so  that  an  equation  of  the 

.7 


form 


-  JL. 

Yt 


<L  <L  <L- 

j  JXj  }  JAj 

the  quantity  6tyK 

Q 


will  be  taken  to  mean  the  vector  has  components 
in  three  dimensions.  Furthermore  by  introducing 
defined  as  follows: 


ev*  =  * 


(1.10) 


if 


are  not  all  distinct. 


-  / 


*//c 


When  the  order  of  cj  / j  K  •  i8  the  result  of 
on  i-j  3  ,  and 


an  even  permutation 


When  the  order  of  Lj  f  j  *  is  the  result  of  an  odd  permutation 

cn  t}  ^  3  ,  the  cross  product  of  two  vectors  can  be  written  in 

compressed  form.  For  example t 


(*  *  ~)i  -  *y**/g«  (1.11) 

Where  the  summation  convection  has  been  used  for  the  indices  J  and  A 
As  a  check,  in  three  dimensions,  the  X  -component  of  A  X  3  can  be 
expanded: 

\  -  €"«  8k  +  e,a/f  bk  /.  tf3  K  bk 
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In  the  first  sum  all  terms  will  be  zero  since  already  has 

repeated  indices;  in  the  second  sum  the  only  non-zero  term  will  occur 
when  X=3,  i.e.  £/lSJt  ^  83  i  and  in  the  third  sum  the  only 

non-zero  term  will  occur  when  K=2t  i«e.  G/m 
Therefore 

(AX  B)x  =  (AgBo-  (1.12) 

This  can  be  seen  to  have  the  typical  form  of  the  components  of  a 


-  _  /  ( *  *  6.  )t- 


Then 


vx£ 


(1.13) 


(l.H) 


(1.13) 


(1.16) 
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Now,  using  the  relations  between  the  fields  given  "in  vacuo*,  equation 
(1,15)  becomes: 


£  *£  s  (1.17) 

taking  the  cross-product  with  j<  , 

£ *(**£)?  *c  xif  (1.118) 

Using  the  vector  identity  Ax(BxC)  -  (a  .  c  -(A-B)  <■  on  tbe  lejft 
and  equation  (1.16)  on  ^he  right, 


£)*  -(«■*)£=  _ 

(1.19) 

from 

-  tu>~r  =■  J  +  ~  ~  J  +  /sjP 

(1.20) 

j)e  =  g  2> 

—  t  \) 

(1.21) 

-  _  €  £  r- 

(1.22) 

(1.23) 

h  --  #. 

where  JT  is  the  identify  matrix. 

In  analogy  with  dielectric  materials,  this  equation  is  written  as 


£e  ^  f.  £<  E  (1.24) 

where  k,  is  called  the  effective  dielectric  tensor  of  the  plasss^ 
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Equation  (1*3)  then  becomes: 


(r'H*  ~  (*  i<)§  =  6  6  § 


(1.25) 


This  is  the  wave  equation  used  in  the  calculation  of  the  complex  refrac¬ 
tive  index.  What  is  needed  for  this  calculation  is  the  assumption  of  a. 
form  for  £  »  the  wave  vector,  and  the  calculation  of  the  dielectric 
tensor  6  •  This  is  where  the  equations  of  motion  play  a  part.  It 

3s  at  this  point  that  one  is  really  justified  in  saying  that  a  generaliza¬ 
tion  of  the  dielectric  tensor  (or  conductivity)  implies  a  generalization 
of  the  Appieton-Hartree  formula# 


B  -  THE  EQUATIONS  OF  MOTION 


The  equations  of  motion  as  already  stated,  are: 


4-  (>»<  Vk  )  =  (£”  +  v*  x  B) 


(1.10) 


but  these,  however,  do  not  take  into  account  forces  which  may  be  introduced 
by  eollisional  processes,  if  or  example,  neglecting  the  magnetic  field  for 
simplicity,  for  an  electron,  (lo10)  becomes:; 


m  1/  -  <  ta  £ 


(1.26) 


where  F0  ,  may  depend  on  position.  This  equation  is  easily  integrated 


i/'  ~  -  (  ( — —  j  €  -f-  6 


(1.27) 
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This  velocity,  it  must  be  remembered,  is  imposed  by  the  electric  field 
and  must  be  distinguished  from  random  velocity*  After  a  collision 


this  'ordered"  velocity  impressed  by  the  field  must  be  zero,  so  at 
time  t/  ,  when  the  last  collision  took  place,  y  z  o  and 


(1.28) 


therefore. 


W  -  /  i 


(B)( 


c'ni'Xi  iu)~C  J 

c  -  e  / 


(1.29) 


If  i  is  set  equal  to  i,  -tl  T  ,  equation  (1*29)  can  be  written  as 


(1.30) 


In  a  senje  equation  (I.30)  gives  the  ordered  velocity  of  an  electron 
in  terms  of  the  time  of  its  last  collision.  If  the  number  of 
collisions  in  a  time  <J(  were  known  the  average  velocity  of  the  electron 
could  be  calculated.  For  ease  in  future  calculations,  an  equation  of 
motion  could  be  constructed  such  that  the  calculated  average  velocity 
of  the  electron  would  be  its  solution.  This  equation  would  include  an 
expression  for  a  "force"  due  to  collisions.  The  times  between  collisions 
are  supposed  to  have  a  Poisson  distribution  f  ,  i.e* 


-ry  / 

/  =  Nyc  (Jr 


(1.31) 
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where  V  is  the  collision  frequency,  then  the  average  ordered  velocity 
is 


e'ICJX)  e~z^ c/r 


(1.32) 


-  <  eVE. 

hi  U> 


AcO-,»)r 

€  1  Sr 
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r  _ /_  7 

yr>u>  1  L*  *»»] 


r  -  U  »£o  aot  /'“EJjUlZ  7 

y»U)  e  L  »(*+*}] 


^ir)  z  tEo  tu/t /_L_  j 

yn  t  (  iV-tp  / 


(1.34) 


If  the  effect  of  collisions  was  to  be  accounted  for  with  a  force  -J(£ 
where  g  is  a  constant,  then  equation  (1*26)  would  be  transformed 
into 


Yn  if' 


(1.35) 


Its  solution  is: 


IT  - 


(tu>fn  +J) 


(1.36) 


If  this  solution  is  to  represent  the  average  velocity  of  the  electron 
as  explained  before,  comparison  between  equation  (I.36)  and  equation 
(1,34)  leads  to 
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-  Irj  y 


(1.37) 


so  that  the  fictional  foce  produced  by  the  collisions  is  • 

The  analysis  here  has  been  simplified  but  a  more  detailed  analysis 
will  change  the  term  by  at  most  a  multiplicative  constant. 

Written  in  full  then  the  equations  of  motion  are: 


4  l-v  Yj  I  -Kc*'?.  (U38) 


C  -  SOLUTION  OF  THE  EQUATIONS  OF  MOTION 

The  subscripts  will  be  dropped.  It  is  understood  that  there  are 

N  similar  equations;  one  for  each  electron.  Furthermore,  all  quantities 

cco  7^ 

are  considered  to  be  time-harmonic,  i.e.  they  contain  the  factor  <?  • 

Therefore,  equation  (1*38)  becomes* 


-  -  -  aom  tfr  +  t  (c  +  cco  r  x  B  ) 


(1.38) 


The  magnetic  field  B  as  previously  stated  was  composed  of  the 
steady  external  field  and  that  of  the  wave  but  it  can  be  shown  that 
tha  magnetic  field  of  the  wave  has  a  very  small  effect  on  the  electrons 
(Ref  12*28)  (See  Appendix  a),  therefore,  it  is  neglected. 
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Setting  3  r  Be  (the  external  steady  magnetic  field)t  equation 
(1.38)  becomes* 


_  f  -  -  1  i*>  >*1  y'  f  -f-  e.  (£  J  c  C  *  -  o  )  (1#39) 

Soffit  n  ^plification  is  possib3.e  if  an  appropriate  coordinate  system 
is  chosen.  The  one  chosen  here  has  its  bz,  axis  parallel  to  Be 
(see  Fig  1).  After  some  rearrangement  equation  (1.39)  becomes: 


te 

u) 


%.*£) 


(1.40) 


where 


,  b 

Us 


because 


>1  «J 


(& o  xrj  -  A..  r 

'A  e*J*  6J  H 


(1.36) 


and  since  B3  f  0  while  .T,  -  -  0  %  as  can  be  seen  by  referring 

to  Figure  (1), 


e  _ 

m  (a) 


<</<  *;  0, 


*<3X  ^3  r* 


(1.41) 


The  ri^it  hand  side  can  be  thought  of  as  a  matrix  multiplication 
between  a  matrix  whose  elements  are  given  by  ^i3K  and  a  column  vector 

fx  7 

(r/t )  -  iij  •  Using  the  definition  of  *  the  elements  of  this 

martrix  can  be  ascertained  to  bs 
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Therefore  equation  (I.38)  can  be  written  in  matrix  form  as  follows:- 


u) 


(1.42) 


At  this  point  several  definitions  are  useful  to  make  some  economy  in 

c  &<> 

writing.  First  ~  s  ,  the  gyrofrequency  of  the  electron,  and 

secondly  rf  -  72  •  Furthermore, 


■'ll: 

yn  u) 


(1.43) 


with  these  definitions  in  mind,  equation  (1.42)  becomes: 


_ i.  [  B„  x  r]  :  -  Jt  t  (1.44) 

>-n 


If  equation  (1.44)  is  used,  equation  (1.40)  can  be  compressed  furthers 


^ (pi*  (A)r  £■ 

0  *  — 


(1.45) 


where  7  is  the  identity  matrix. 

It  is  understood  that  this  equation  is  a  matrix  equation.  It  can  be 
seen  easily  that  [pX  /  ;A  J  i3  non-singular,  therefore,  it  has  an 
inverse,  and  the  solution  to  the  equations  of  motion  is 


(1.46) 
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D  -  THE  DIELECTRIC  TENSOR 

Tlie  solution  to  the  equations  of  motion  just  obtained  is  not 

very  illustrative  in  itself.  Its  main  usefulness  lies  in  the  fact 

that  from  them  the  conductivity  and  dielectric  tensor  can  be  calculated* 
It  tvill  be  recalled  that  this  is  the  problem  which  the  Lorenta  theory 
promised  to  solve*  and  in  turn  the  determination  of  these  constitutive 
parameters  makes  it  possible  to  determine  the  index  of  refraction  of 
the  medium  from  the  wave  equation  (1.25)* 

It  will  be  recalled,  that 

J  =  6  £  (1,7) 

and  also,  that 

J  -  P  (1.47) 
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WVta 


?  2 

, }  „  3I--  Is  the  square  of  the  piasna  frequency*  It  is 

*  "  t»6* 


cleai1  frco  equation  (1.49/  that  the  conductivity  tensor  is  given 
by 

i-/ 


6'  - 


A)' 


(1-50) 


4 

% 

4 

* 


Since  the  dielectric  tensor  is  defined  by 


■hi 


(1.51) 


it  is  clear  then  that 


cc  r  ~  &(/$!  +  ! A)  h  1 


(1.52) 


where 


When 


C/3f 


z  s  J  +  T 
v 


-?  j  - 1 


is  set  equal  to  ^  ,  then 


U.53) 


3  is  called  the  susceptibility  matrix. 

At  this  stage  one  could  proceed  to  calculate  a  dispersion  relation 
from  the  wave  equation  (1.25)  aa&&  therefore  get  the  index  of  refraction  of 
the  medium.  It  is,  however,  preferable  to  have  in  diagonal  i'orai* 

that  would  certainly  simplify  the  dispersion  relation,  which,  as  a  rule, 
is  always  fairly  complicated.  This  latter  procedure  is  done  by  every 
author  consulted.  Briefly:,  it  consists  of  finding  &  non-singular  matrix 
p  such  that 


A  At  A  = 


(1.54) 


?3 
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where  is  diagonal.  Thin  is  called  a  similarity  transformation 

-  w  —  ' 

•y 

and.  is  said  to  be  similar  to  €e  •  From  equation  (1.53)  it  can 
be  seen  tbst  the  same  matrix  A  which  diagonalizes  »  also 

diagonalizes  S  and  hence  C  since  -5" -  «  3y  the  same  token 

if  it  diagonalizes  g  it  diagonalizes  £3/  +  i  Ji  )  .  And  since 


(pl-f  iA)  1  (px+t  A)  ~  T 
Af(pI-f(Aj  'AA'* (fil+iA)  -  A  '/ 


(1.55) 


(1.56) 


D  A  '‘(pi A  (A)*  -  1A ' '  (i.57) 

inhere  d  is  a  diagonal  matrix  by  hypothesis.  If  both  sides  of  equation 
(1.57)  are  multiplied  by  A  from  the  right,  equation  (1*57)  becomes 


DA  'f  (pj  -hA)A  -  /  (1»58) 

If  it  is  further  supposed  that  the  elements  in  the  diagonal  of  i>  are 
all  non-zero,  then  J)  is  invertible,  and  hence 

fi-'(p T-finJ/i  -  r>"  <1-595 

But  the  inverse  of  a  diagonal  matrix  is  also  diagonal,  hence  it  has 
been  proved  that  A  also  diagonalizes  (jslA  1  •&  /  • 

There  is  a  theorem  in  matrix  algebra  which  asserts  that  if  the 
eigenvalues  of  a  matrix  are  all  distinct,  then  the  diagonal  matrix 
with  the  eigenvalues  along  the  main  diagonal  is  similar  to  the  given 
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-nv  /' 


>1 


matrix.  Hence  if  the  eigenvalues  of  (  fil+tJi  )  are  determined, 
and  if  they  turn  out  to  be  distinct,  D"^  can  be  written  down 
automatically. 

y\ 

The  characteristic  equation  of  (  jbZ  +  )  is:. 

-iA 

/  *  0 


if  '  A)  l  (/3-j)~~  SlL]  =  0 

from  which  the  eigenvalues  are  seen  to  be: 
A,  sr  P+-A- 

K  -  p-i*- 

\3-f> 

Therefore , 


since  all  the  eigenvalues  are  distinct. 
From  equation  (1*59) 

(pI-hiA)A  =  At)'1 
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(1.62) 

(1.63) 

(1.64) 


(1.65) 


(1.66) 


i 

i 


I 


! 

i 


Not  all  equations  in  this  group  are  independent,  so  that  not  all  the 
elements  of  A  can  be  determined.  This  merely  means  that  the  matrix 
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jC  *  * +  */  v  ^  ^ 


A,  is  not  unique.  The  relations  that  can  be  obtained  are  the 
following : 


=  +1 


*31  =  ° 


a12  - 
^22 


=  Arbitrary 


(1.68) 


(a  +  bi)(a  -  bi)'  =  a2  +  b2 


(1.69) 


The  elements  in  columns  of  the  matrix  A  are  frequently  interpreted 
to  be  eigenvectors.  Their  magnitude  should  be  1.  Therefore: 


(All|2  +  fail  2  +  | A3i (  2  =  l  (1.70) 


|/in  L  i 

r  il Zi* 


(1.71) 
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Z  = 


M// 
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(1.72) 


KJ'-  7= 


/r 


(1.73) 


And  also 


IA„  I  -  j= 


(1.74) 


Similar  equations  hold  for  the  second  column,  and  for  the  third 
A 1,  Thus,  as  a  matter  of  convenience,  the  matrix  here  chosen  is 


*r 


111  o  i 

-i  i  0 

V  0  0  /2 


(1.75) 


The  inverse  of  this  matrix  A  car,  be  calculated  in  a  straightforward 
manner  and  turns  out  to  be:; 


A' 


-1 


1 

n 


I  1  i  0 
1  -i  0 

Vo  0  J2 1 


(1.76) 


It  will  be  observed  that  this  matrix  is  also  the  conjugate  transpose 

*1 

of  A  denoted  by  A  .  In  other  words,  the  matrix  A  is  unitary. 
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The  intent  of  these  preliminaries  was  to  obtain  A  so  that  the 
diagonalization  of  could  be  effected,  and  thus  simplify  the 

final  dispersion  relation. 

It  is  recalled  that 
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/'  0  *  I 

\  '  0 
{  6  6  /  / 


(1.78) 


*  -  77= 


(1.79) 


1  *  ^  .  «  / 

0  I-J, 

To  reduce  the  writing,  the  following  definitions  are  made:- 


x-  &- 

'  ft + 


/3  4-JTL 


(1.80) 


-  /  - 


ft- n 


(1.81) 


^  =  /- 


(1.82) 


so  that 


= 


(1.83) 


Finally,  to  transform  the  ware  equation  to  this  new  system,  where 
fc't  is  diagonal,  the  transformation  law  for  the  vectors  must  he 
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obtained.  This  may  be  obtained  .from  equation  (l.?4). 


4  -  €0  ic  0 


(1.24) 


If  the  equation  is  multiplied  by  A“  from  the  left,  then 

A”K  =  §  =  *.(*%/! J*'f£  = 


(1.84) 


'l  =  t.UE 


(1.85) 


where  De,  E  denote  De ,  E,  in  the  new  system;  therefore  the  transforma¬ 


tion  ldw  for  a  vector  X 


=  A'1! 


(1.86) 


Vfaere  the  bar  on  top  denotes  the  vector  in  the  new  system.  In  the 
present  calculations,  equation (1.86)  represents  the  following  three 


equations: 


khh  +  i'lhl 

®2  ~  ^[^1  "  1  (—^2 1 


(1.87) 


(1.88) 


(a)^  =  (X)^ 


(1.89) 


These  transformation  equations  are  of  great  importance  in  the 


calculation  of  fields  in  a  plasma. 
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The  final  calculations  on  the  wave  equation  can  now  be  made. 
If  in  the  wave  equation  (1.25)  the  wave  vector  k  is  taken  as 


(1.90) 


where  /  is  the  speed  of  the  wave  and  M  is  a  unit  vector  in  the 


direction  of  propagation,  the  equation  becomes 

(7!  (s-S)s- 


(1-91) 


(P'i)*  ~  (-'-)£  -  (1.92) 


If  this  equation  is  transformed  by  4 


(to  -  §  )  Yi  _  S  =  -  Me  V*l>e 


.  (1.93) 


where  it  has  been  kept  in  mind  that  the  inner  product  of  two  vectors 
is  invariant  under  linear,  homogeneous  transformations.  The  second 
term  on  the  left  of  (1.93)  can  be  written. 


r-i  - 

( t );  '  —— 

(  ("e  )( 

where  L  denotes  a  particular  component  of  the  vectors, 
Therefore,  in  component  form  (1.93)  can  b e  written  as 


(1.94) 


(yj-t) 


(1.9 5) 
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In  general  ( fr  ■  t )  is  not  zero*  thus  avoiding  the  trivial 

v 

results  which  would  otherwise  be  obtained* 

If  the  inner  product  of  the  terms  in  equation  (l.l6)  with  n  is 
taken  (i.e.  with  <  )  then:. 


-  K  •  (*  *  H 


(1.96) 


Or  (J 


K  .  A, 


(1.97) 


So  that  in  this  case  JJ  is  orthogonal  to  J>(  -  Assuming  that  2?  is 

real,  then  n  is  also  orthogonal,  and  finally  because  a  unitary 
transformation  preserves  the  complex  inner  product,  it  is  also  true  that: 


h*.  Dc  h*.  4 


e  --  O 


(1.98) 


where  n  *  is  the  transformed  vector* 

Iherefore  multiplying  (1.95)  by  n  *  ,  the  result  is: 


—  -  K2 


— 

^  4  * 


-A- 


(1.99) 


or  using  equations  (1.87),  (1,88),  (I.89): 


~  M.-,  V 


r7  -  A*.  V ' 


(1.100) 


r$$pts>»V ’aw*  V’ 
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t(”Mj 

-x-  .  1 /> 


O 


(1.101) 


and  defining 


f2 


(1.102) 

(1.103) 


/3  *-  —7 —  (1.104) 

*<>  fj 


the  equation  becomes* 


(1.105) 


luis  id  the  dispersion  equation  which  gives  a  quadratic  aquation  for 
yl .  This  means  that  under  the  conditions  assumed  in  this  derivation 
there  are  in  general  two  waves  which  propagate  at  different  velocities. 
These  waves  are  called  the  ordinary  and  the  extraordinary  depending  on 
how  much  they  are  affected  by  the  magnetic  field.  Thuse  two  modes  of 
propagation  are  due  to  the  effect  of  magnetic  field.  It  will  be  seen 
that  in  the  final  formula  if  the  magnetic  field  is  set  equal  to  zero, 
any  distinctions  between  the  two  waves  vanish. 
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D  -  THE  AFPLETON-H.-^PT'R£E  EQUATION 

Supposing  now  that  the  direction  of  propagation  is  entirely 
arbitrary  (Fig  1)  then  ^  can  be  written  as 


H  -  (sin $  cus  Si)i&s/r><tt  cos  &  J 


(1.106) 


in  a  spherical  coordinate  system;  one  can  see  that  ft,  V si”2'  & 
and  ^3  '  C^s  &  so  that  it  is  evident  from  (1.105)  that  the  equation 
resulting  from  (1.105)  is  independent  of  £  ,  but  not  of  6  ,  which*  if  it 
is  recalled  that  the  direction  was  chosen  parallel  to  the  magnetic 
field,  is  the  angle  between  the  direction  of  propagation  and  the  magnetic 


field. 

Equation  (I.I05)  can,  therefore,  be  written  as: 


J-  (s  '/>j  ) 

^7r 


cos  X& 


(1.107) 


Furthermore,  tne  complex  index  of  refraction  is  equal  to  -~- 

where  \fc  is  some  reference  velocity  like  the  velocity  of  light,  so  that 


(1.108) 


and,  equation  (1.107)  can  be  written  as: 


2. 


&/*)*& 


2 


cos  '^6 


0 


(1.109) 
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Now  since 


e- 


(1.110) 


o  o 

and  if  Vo  -  c^  equation  (-109)  can  then  be  written 


s;»2e  .  HrH  +■  lc0,s*e 


6,  j-l1 


i _ 


i _ L 

6*  yiv 


(1.111) 


Prom  this  equation  the  Appleton-Hartree  expression  for  the  complex 
index  of  refraction  can  be  obtained.  (See  Appendix  S).  The  result 


2  i  -  v  a 

M  r  i  ~  ,-tZ-  —L. 


AI - -  +  /  >r _ . 

iO-X-iZ)  -  J4(i-X-iZ)x 


v: _ „ 


(1.112' 


where 


KrO. 

» -iZ-7 
Yr  A 

Yr  r  Ks/yifi 

^  r  ^5  0 


and  can  be  written  as  (^  -  //  )  =  R.~cT  .  The  positive  tarn 

denotes  the  ordinary  wave  and  the  minus  sign  the  extraordinary  wave. 


3P/PH/C5-1 


It  is  easy  to  see  that 


M  '  /  R  +  1/ z2 


(1.113) 


f" 


(1.114) 


for  each  sign.  Hence  it  follows  that  both  ^  and  ^  are  always 
non-negative,  and  from  the  fact  that  the  field  quantities  contain 
i(ej t- £.'£)  iUitt  ~  ~  r  ) 

the  factor  e  ~  £  where  it  ha3  been  assumed 

that  C  ,  then  substituting  for  M 

r  ■ Wr]  t&t- 

e  =  &  (1.115) 


Mt-gjur)  -*Kr 

v-  a 


(1.1-  *) 


as  it  is  seen  that  apart  from  an  oscillating  factor  C 


U*>t-T*r) 


all  field  comnonents  contain  a  factor  £  c 


rf  r 


,  which,  because 


is  positive,  decreases  as  r  increases.  This  lerm,  therefore,  represents 
an  attenuation  on  the  wave  as  it  travels  in  the  direction  of  increasing  r, 

~  X  is  defined  as  the  absorption  coefficient  of  the  wave.  The 
preceding  analysis  has  been  concerned  with  only  a  point  in  wbat  is 


'  '  -  •%**- 
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presumed  to  be  a  homogeneous  medium.  In  the  ionosphere,  however, 
both  the  electron  density  and  the  collision  frequency  are  functions 
of  height  (and  if  one  wishes  to  be  more  exact,  so  is  the  magnetic 
field),  so  that  the  expression  given  here  for  the  fields  are  not 
strictly  valid,  for  account  must  be  taken  of  the  variations  mentioned* 
The  ionosphere  is  visualized,  in  this  new  approach,  as  horizontally 
stratified,  i.e*  N  and  y  do  not  depend  on  if  Z  has  been 

chosen  as  the  height.  This  means  that  M  is  a  function  of  height,  so 
that  the  differentiations  made  at  the  beginning  must  take  into  account 
this  fact*  This  is  done  by  Sudden  (Ref:  12),  The  solutions  are  the 
so-called  WKP  solutions  of  the  ionosphere.  For  the  particular  case 


of  a  wave  horizontally  polarized,  i.e. 
component  is  t 


the  solution  for  this 


(1.117) 


for  an  upgoing  ware.  Upon  substitution  of  Ms  /*-*)( 


(1.118) 


and  assuming  that  the  variation  of  M  wita  height  is  fairly  small,  it 
can  be  concluded  that  the  factor  which  indicates  the  attenuation  of  the 


wave  is 

c 


(1.119) 
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or,  in  terms  of  the  absorption  coefficient, 


-/? 


(1,120) 


The  integral  is  the  absorption  cf  the  wave.  It  is  this 

integral  v/hich  is  the  subject  of  calculation  in  this  paper. 

Now,  since  the  Appleton-Hartree  equation  refers  to  a  single 

point,  no  mention  being  made  of  how  M  varies  from  point  to  point, 

it  is  then  assumed  that  it  is  valid  at  any  point.  At  each  point, 

then,  the  proper  values  for  the  parameters  are  substituted  in  the 

2 

equation  to  obtain  the  proper  value  of  M  .  In  principle  f 
can  be  calculated  then. 

The  condition  for  validity  of  the  WKB  solutions  has  also  been 
given  by  Budden.  This  condition  is: 


_/  J—  j  <<  / 

/ q  [sf2-  JIJ  Zz~r 


(1.121) 


It  is  evident  that  this  condition  fails  when  tf-0  ,  i.e.  near 


levels  ef  reflection.  The  derivative 


c/lM 

and  should  be 


fairly  small  for  the  condition  to  be  valid;  given  the  term 
the  condition  is  easier  to  obtain  at  higher  frequencies* 


sp/ph/65-i 


Gexieraiitt&iioiis  uf 


xi.  rts  _ _ j  .1  mu  _ 

mi e  incur)' 


A  more  general  approach  may  be  taken  to  calculate  the  index  of 
refraction  of  a  magneto-ionic  medium.  The  most  iramediate  step  in 
generalizing  the  theory  presented  in  the  previous  Chapter  is  to 
consider  more  carefully  the  collision  of  the  electrons  with  the  other 
constituents  of  the  medium.  It  was  implicity  assumed  in  the  classical 
applston-Hartree  theory  that  the  collision  frequency  of  the  electrons 
with  the  neutral  particles  of  the  magneto-ionic  medium  was  constant. 
Experimental  results  have  shown,  however,  that  for  gases  in  general 
this  is  not  true.  The  electron  collision  frequency  turns  out  to  be 
velocity  dependent.  In  particular,  in  nitrogen,  it  is  proportional, 
to  the  square  of  the  electron  velocity. 

The  Sen-Wyller  theory  then  is  the  generalization  of  the  Appleton- 
Hartree  equation  to  the  case  of  non-constant  collision  frequency.  The 
techniques  used  are  similar  with  one  major  exception:  in  the  Sen-Wyller 
generalization,  it  becomes  too  difficult  to  deal  with  the  equations  of 
motion  directly  so  that  use  is  ma*.e  of  Boltzmann’s  equation,  which 
permits  the  calculation  of  number  densities  and  average  speeds  once 
the  electron  distribution  function  has  been  found. 


A.  The  Assumptions 

In  the  first  place  the  medium  is  assumed  to  be  a  slightly  ionized 
Lorentz  gas.  This  Lorentz  gas  two  constituents  of  more  importance 
than  others  that  might  be  present.  These  are  the  electrons  and  the 
neutral  particles.  Furthermore,  the  mass  m,  of  the  neutral  particles 
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is  much  greater  than  that  of  the  electrons  .  The  slight  ionizsLtioa 
condition  cen  be  expressed  by  where  H,  ,  is  the  number 

of  the  neutral  particles  and  yiv  that  of  the  electrons.  It  is  also 
true  that  H£  ,  the  collision  frequency  of  the  electrons  against  ions* 
is  much  less  than  F£/l  ,  the  co)lision  frequency  of  the  electrons 
against  neutral  particles.  This  is  what  is  meant  by  saying  that  the 

electrons  and  neutrals  particles  are  the  most  important  constituents. 

-  -  '  . 

of  this  medium.  V  is  supposed  to  apply  to  elastic  collisions. 

The  plasma  is  assumed  to  be  uniform  and  under  the  influence  of 
a  constant  magnetic  field  and  an  oscillatory  electric  field. 

Under  the  conditions  described  above,  several  preliminary 
inferences  can  be  made  which  help  to  solve  Boltzmann's  equation. 

B  ~  Preliminary  Inferences 

f'or  particles  with  a  distribution  function  (the  electrons  in 
the  present  derivation),  the  general  Boltzmann’s  equation  is 


if  *£■  %.  t  =  &  (2.1) 

v  T 

^"here  F  is  the  force  per  unit  mass  on  the  electrons,  v'z  their 
velocity,  Vy  the  gradient  w^tor  in  velocity  space,  ’*>  the  gradient 

X 

in  position  space  and  £  is  the  collision  term.  Under  the  assumption 
of  plasma  uniformity  fj,  should  be  independent  of  position.  Therefore 
a  4  ~  0  ,  nd  equation  (2.1)  simplifies  to 

vr  *  z 

f 

S  /  f  /*| 

^4  +  f  .  K  =  u 


hi 


(2.2) 


lyr'/C-i 


F  can  be  easily  calculated  since  it  represents  the  external  forces 
acting  on  the  electrons  per  unit  aass.  The  m  £n«tic  force  is 

X  Ho  *  ~d  the  electric  force  is  e^Tcosvt ,  ,jh?re  2  is  time- 
independent  .  Therefore 


-  (eoI2X  K_ 

O}  -  L  “N. 


•*  eJSccswt) 


G.  c  oswt  -  v  /.  X 

_  —5  —O 

r**»  _  *■ 


(2.3) 


./here  i  2  = 


—  £ 
vr,  - 


,  Hence,  if  equation  (2,3)  is  substituted  in 


e_ action  (2.?)  then  the  1  -tier  becomes 


4.  (p  cos  lot  -f  _*£  v  X.  H, )  •  Vv  "fi  -  Q 


(2.4) 


The  last  term  which  sho  Id  be  examined  is  the  collision  term. 

This  -.:ill  be  done  only  briefly  sin-e  any  calculation  of  4  is  fairly 
involved.  The  reader,  ini  crested  tn  these  details,  •  ouli  do  well  to 
consult  Chapman  and  Cowlinc  (nef  15043).  In  pener-l,  it  can  be  said 
that  if  there  ~re  two  constituents  an  the  plasma  :itn  distribution 
functions  f^  and  ,  and  if  the  collisions  between  particles  of  these 
constituents  are  elastic  with  relative  velocity  r  at  the  time  of 
collision  then 


Q  -  [{[(// fa  “  iD^dU&Jv', 


(2o  5) 


are 


,  i-  i  Vi  ..  <  t  -j-  'x 


The  primes  indicate  the  distributions  after  collision,  ^  and  f 
impact  parameters  analogous  (though  not  the  same)  as  those  used  in 
Rutherford  scattering.  These  definitions  may  be  found  in  Chapman  and 
Cowling  (Ref  s  15).  Finally  (Jv>  stands  for  (c/lC )x(d^ (dd  )z  , 

that  is  an  eleaent  of  volume  in  a  3“dimansional  velocity  space. 

One  final  inference  and  the  equation  (2.4)  should  be  in  its  most 
appropriate  form.  This  inference  is  that  should  be  a  Maxwfllia® 
distribution  since  the  mass  of  the  neutral  particles  is  much  greater 
tuan  that  of  the  electronc,  i.e. 


1  '■  ».  ( !rTr/ 


I  */*  -  Cl. 


/  X 


'•here  -<  ,  is  tne  number  density  of  the  neutral  particles  and  the  other 
symbols  have  their  usual  meaning.  This  is  the  reason  why  only  one 
Boltzmann's  equation  has  been  written  down.  Normally  a  Boltzmann 
equation  would  have  to  be  solved  for  each  constituent  of  the  plasma. 


C  -  The  Solution  of  Boltzmann's  Equation 

The  method  of  solution  is  the  so  called  Ohapman-Snskog  Method, 
’’heoretical  justification  for  this  method  may  be  found  in  Chapman  and 
- *•  1 1  rg  (.oef  15)  »  Here  as  in  the  Sen-Wyl) er  paper  (Ref  49)  it  will  merely 
re  used.  Tne  method  consists  of  the  following,  Jz  is  assumed  to  be  of  th 


/  V  ;7l  •  14 /(V  cosiJT  +pxS/7teJt )  (if,  xfftj-Y,  C^ccjOi  + 

-  >1  id  x  {d  x3)j'yz  //*  Cc-sc-lt-i  SiSirU*s? ) 


(2.7) 
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r  101 

fl here  7  '  /j  /  /  /,  are  dependent  only  on  the  speed 

)  *)  1  (-  j  1 ; 

of  the  electrons.  The  task  of  solving  Boltzmann’s  equation 
reduces  the»  to  the  determination  of  the  functions  just  mentioned. 

If  equation  (2.7)  is  substituted  in  equation  (2.4)  and  the  coefficients 

of  i\  ^  / \  .  //4  /V,  X  /I  J  V,  and  CtL'  *  (fl0  jl  '  on 

both  sides  are  equated  to  each  other,  then  with  suitable  manipulations 
(see  Appendix  g)  a  system  of  linear  simultaneous  equations  results. 


From  it  0/ 


j,  ¥  m  Y“  J  can  be  calculated.  The  results  are: 
)  r*j  K  j4,  ** 


<*  -  -  ^ 


(2.8) 


P*-  *  V*-/-  ^ 


(2.9) 


_«i  (yl+s*-vt’)  _ _ _  j. .  j/,?* 

^  [» V  (u)*sp]l»\  (cos)*]  %  Jir 


(2.10) 


>1  -  __r 
"  »» 


/ . 


z  [yV(*)+s )*][»*+ (vsj* ]  <5. 


(2.11) 


/  . 


///  />  V  (V*Jp][y  v  (<J'S)X]  ^ 


(2.12) 
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r  -J.1  _/ _  ufs*--  +  )  _ ,  • 1 .  *  (2.13) 

Hf  (»*■**>*)  /iv  ^  ^ 


The  symbols  have  been  defined  in  Chapter  1.  Only  f 


remains 


to  be  calculated.  Th  3  may  be  done  from  the  following  equation; 


(2.14) 


V/here  ^  is  tne  angle  between  l\  and  f£o  and  \  )  is  the 

mean  free  path  of  an  electron  with  speed  ^  .  The  mean  free  path 

can  be  expressed  in  terms  or  collision  frequency  as  * 

If  equations  (2.8),  (2.12)  are  substituted  for  t  and  K. 

in.  (2.14)  then 

nz  3  (fsJi-  i-  X  _ _  _i_ .  ;£o) 

<x  Jirx  J  V*4ut£v^(uj+s)2][vt+f0!  !>)1]  ^  <?*4 


nr5fiiW+ ^V’ 

ri,  X  x  ^  yvj,  A  2 


(2.15) 


If  C  is  defined  to  b' 


52-'v)>  4  6*-  3U>*  ) 

[^Oo+sf]^*  fa-*)*] 


(2.16) 
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and  if  a  “  -j  is  used,  then  equation  (2.15)  become* 


**».  P14 


(2.17) 


or 


rf(K)' 
t  (vN/). 


-  K1^ 


k' 


0) 


(2.18) 


and  if 


t  (y*+  u1) 


then 


ii('K 

—  (  KT  +  F(*i)) 


4  u;  -  0 


(2.19) 


This  differential  equation  which  is  linear  with  non-constant 
coefficients  can  be  easily  solved,  i.e 


KT  + FOrJ  J 


(2.20) 


B  is  jast  the  normalization  constant.  Several  interesting  facts 
emerge  from  equation  (2.20).  First  of  all,  not  even  f_^  is 

r, 

Maxwellian.  The  deviation  from  a  Maxwellian  distribution  is 
controlled  by  size  of  F^g).  This  equation,  it  should  be  stressed. 
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is  valid  for  any  slightly  ionized  Lorentz  gas.  As  will  be  seen,  this 

/  c») 

will  be  specialized  later  to  the  ionosphere  where  /,  will  be 

assumed  to  be  Maxwellian  which  implies  certain  relationship  between 
the  size  of  r(^x)  and  K T  . 

Without  writing  j ^  in  full,  it  can  be  said  that  Boltzmann' 3 

equation  has  already  been  solved  since  the  functions  ^  /J,  jr*  h  ^  £ 

r(t)  * 

and  A  of  ^  are  dow  known. 

From  here  on  the  procedure  is  largely  similar  to  that  in  Chapter 
I.  The  next  step  is  to  calculate  the  conductivity  and  the  dielectric 
tensor . 

D  -  The  Conductivity  Tensor  and  The  Dielectric  Tensor 

It  may  be  recalled  that  in  the  previous  chapter  the  conductivity 
tensor  was  obtained  by  deriving  an  explicit  relationship  for  f ■=  € t  . 
This  relationship  was  based  on  the  solution  to  the  equations  of  motion 
jr  i  /"  (£ )  .  In  the  Ser-V/ylier  generalization,  no  such  relation 

sxi  st3  but,  however,  since  /  where  <¥.*)  is  the  average 

velocity  of  the  electrons,  and  sines  (V  sad  <J4)>  can  be 
calculated  f~os  the  knexn  distribution  function  4  ,  then  J.  can 
be  calculated  and  presumably  the  calculated  J  can  be  set  equal  to 
£  t  and  o  calculated. 

For  convenience  is  written  as: 


/J'1  *  //l 


(2. 21) 
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so  that  one 


/l  fc  cosidt-f- px  shicol)  +  (H*  x/1  )(&  usvt + jf> 

+  [H,  x(M*  Y-{\  to$u)t+  <JT  si»u)t ) 


5,y:  toT) 

(2,22) 


The  number  density  of  the  electrons  is 


M 


4  -  /  A  J* 

4  LS 


(2.23) 


‘.■/here  J\J  is  an  element  of  volume  in  velocity  space.  Since  ^ 

depends  on  the  speeds  J~2  it  is  used  to  express  JV  in  a 

spherical  coordinate  system  in  velocity  space.  Hence  </?-  ^sihddsajt/di 

Again  for  convenience,  the  spherical  system  here  chosen  is  a  special 

jOl 

in  that  it  is  chosen  such  that  tx  is  parallel  to  the  9  -  0  direction 
(see  Fig  l.b).  Therefore: 


V  =  J(K0)4-  f2°-Vx  )^s/»  e  </*  </j  ^ 


(2.24) 


or 


^  iff  ,tf  r*C 

Nx'  4nf  -h  f  cty  °sin* cos 6 o'6  / it3 fa 

O  *  Jn  e>  "a  *  * 


(2.25) 
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so  that*  finally* 


/£  4!ij 


(2.26) 


Interestingly  enough  then  /%,  is  determined  strictly  by  . 

This  equation  also  gives  the  normalization  condition  for  the  constant 
3.  Using  the  unnoi'malized  distribution  function  then* 


K  >,  =  /<?,  L 


(2.27) 


but,  in  the  spherical  system* 


sib  6  u>S  d 


(2.28) 


so  that 


Jussive  cos  $  )(£  (c)f  cos  6)(/3sib60,6J4JiTtL  (2.29) 


2// 


Both  integrals  are  o  >  since  /  cos  4  ^4  -  <2  .  The  situation 

is  similar  for  the  y-coraponent .  For  the  z-component,  however. 


~  jtK  COS  &  (fz  }+  £  cose  loosing  (f $  <ff  (fc,^  (2.30) 


-A  tc>oih;„acos  d  </ecf4<f<r2  ^ 


(2.31) 
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The  first  integral  is  O  since  f  C6S&  c/e  -  O  but 

the  second  one  is  different;  in  fa,’.'  ,  i.  equal  to  j ~  f cf <^2. 

It  is  clear  then  that  J  is  parallel  to  ^  »  so  that  if  the 
spherical  system  used  were  a  move  general  one  in  which  had  an 

arbitrary  direction  then 


(2.32) 


and*  therefore 


I  = 


47/ e2 


f  fj\  V<£ 


C 2.33) 


3  4 


so  that  -4  determines  S  .  The  partial  objective  stated  at  the 
beginning  of  this  section  has  been  accomplished;  that  is  to  calculate 

*£?< /7/VA. 

o  A  — * 


3  o 


y  .  Now,  the  final  step  is  then  to  express 
as  6  €  . 

If  equation  (2.22)  is  used,  equation  (2.33)  becomes: 

/-? Jj«,  c**t*/*, 

+  [&<,  *(&*!})]/'  (r%0>sa?+S,s/>futJtr/cji 

0  (2.34) 


or 

~  r  'iJ'£  *  Jx  (#•  x£*)+J3  l\  x(K  x£*j] /  (2.35) 
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where 


Jt  r  /  (°<i  uscjf  •/ px  s/'wtjcr 


(2.36) 


/'c°  / 

4  r  4  ($*  cos  tot  +  s/»a)t)or  *  Jo \ 


(2. 37) 


j3  -z  j  (Px  cos  tut  -/-  $x  s 'r>  tot)  ^ <J(Fl 


(2.38) 


A  coordinate  system  must  be  chosen  to  calculate  6  »  The  coordinate 

system  chosen  here  is  similar  to  the  one  in  Figure  1  with  the  exception 
that  the  propagation  vector  ?  is  in  the  x-z  plane  (Fig  l.c) 

If  use  is  made  of  the  notation  developed  in  Chapter  I,  then 


(  ll°  *-)ij  S 


From  Figure  l.c  Hx  =  -  0  }  M*,  =  Ho 


(2.39) 


,  so  that 


=  iit  kMC-4 


(-.40) 


fis  was  done  in  Chapter  I,  equation  (2.40)  can  be  interpreted  as  the 
product  of  matrix  (  H3  )  ctnd  a  column  vector  -  —  (-)*  * 
Therefore,  by  the  definition  of  tjsK.  ♦  in  matrix  notation  equation 
(2.40)  becomes 


-  W.  n  1/  f.  J 


/  0  n-  ? 

°[e/] 
1  '  0  0  0  /  \  t}  / 


(2.41) 
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it  is  clear  that 


-W3  0  |/(H0x£\j|  I 

0  O  //(H.xrj  I 
0  °'WM.xriJJ/ 


and  by  (2.4l) 


0  "  Hj 

0  1/  O 

-M3  oj 

M***)] -Sr.fr  0 

°b 

0 

°  I 

0  yi  0 

£ 

0  /\ 

(2.42) 


(2.43) 


(2.4/0 


so  that  by  equation  (2*35)»  (2.4l)»  (2.4 4}  then 


3 


J,H, 


o 


-  j»h3 

0 


o  |  E*  /  (2.45) 

j,  /( 


30  that,  finally 


-  *fll  '-l 


'I 

•M. 


3  »>. 


-JIM. 

J,  -Me*J3 

o 


.  /  X 


i 
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The  dielectric  tensor  can  be  calculated  easily  enough  from  the 
*• 

definition  £  z  Z  +  t—t  «  To  simplify  the  writing  somewhat 
the  following  definitions  are  useful: 


JjL 


/ 


and 


(Tr  ~ 


/- 


4ffi  u>Z  j' / 
3  CO  ~ 


(2.47) 


3  co 


(2.48) 


*  CO  * 


/ 


Using  the  previous  definitions,  then  fe  becomes 


(2.49) 


o 


(2.50) 


This  was  the  main  objective  stated  at  the  beginning  of  this  section. 

A 

The  comments  that  appeared  on  Chapter  I  on  the  use  of  €c  as  it  stands 
at  the  moment  in  the  dispersion  relation  also  apply  here.  In  Chapter 
I  it  was  diagonaliaed  and  then  the  dispersion  relation  was  obtained 
using  the  diagonal  •  To  show  that  that  is  not  the  only  simplified 
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procedure*  in  this  chapter  a  less  general  one  is  used  but  which  gives 
the  same  results. 


E  -  The  Dispersion  Pelatioa 

It  will  be  x’ecalled  that  the  wave  equation  is 


A 

(n  *f  )J?  ~  (n  *n)£  s 


(2.42) 


where  w  is  a  unit  vector  in  the  direction  of  propagation.  If  the 
propagation  is  parallel  to  one  of  the  coordinate  axes*  then  (f?  # 

will  certainly  be  simplified.  This  is  the  basis  for  the  procedure 
followed  by  Sen-’.vyller  (Ref  49)  and  this  is  also  the  reason  for 
choosing  the  coordinate  system  in  Figure  l.c,  for  then  if  the  system 
in  Figure  l.c  is  rotated  about  the  y-axis  as  to  bring  the  propagation 
vector  to  a  direction  parallel  to  the  x-axis,  then  the  rotation 
matrix  of  the  transformation  is  considerably  simplified. 

The  rotation  matrix  is  (see  Figure  l): 


A: 


o 

/ 

0 


(1.92) 


where  the  angle  ~  $  where  p  is  the  angle  between  the 

magnetic  field  and  the  propagation  vector.  The  matrix  A  is  real  and 

4  A/ 

orthogonal*  therefore,  the  dielectric  tensor  in  the  new  system  is 


3?/P K/65-l 


where  is  the  transpose  of  A.* 

Therefore 


^  /COJ  4 


s/»4  \/tT4t7r, 


/ 

(0 


o 

COS$f 


*3 


0  o  -st*dj 

*x+£n  6  I  o  ,  o  7(2.52) 

c  €x/[*s;»i  0 


or 


'(?  +  <^7?  Cos1  4 

61J  COS  4 
-  6/7)  si*  4>  cosj 


-  eg  cos  4 
*1+64 

ess'»f 

The  wave  equation  (2.92)  leads  to  three  equations 
—  (*l  6  6/22  ^  4)- £%  -b  6jj  Cos  4  Ej  -f  6/n  3/^4  cos ^  £*.-() 


(2  .53) 


(2.54) 


-  to  usf)£,  --  o  (2_55) 


f:rn  sitjjtcostf  Fx  -  6g  si» <f>  Fj  /  [** (til tir, X'^jjF^  -  0  (2.56) 


where  by f  fz 


are  the  vector  components  of 


£ 


in  the  new  system. 


*  This  seems  to  be  contrary  to  what  was  done  in  Chapter  I,  where 
it  would  seem  that  (t  in  the  new  system  would  be  Jt*9A  .  there  is 
no  contradiction  if  it  is  observed  that  the  term  rotation  matrix 
applies  to  the  vectors,  i.e.  €***&  ,  where  £"'is  the  vector  /fin 

the  new  system  so  that  the  similarity  transformation  on  is 
actually  (W7  )  ’ef  (  A’1  )  or  A(e/f  in  this  case. 
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The  astern  of  equatic  *s  (2.54) »  (2«55)  and  (2.56)  has  a  non-trivial 
solution  only  if  the  determinant  of  the  system  is  0  .  Therefore, 


~(tr+6/W 

6n  CPS<f> 

tg  sin <j>  cost 

6jt 

eg  sintcosj 

-  6jt 

MK  (€t Umsinl4>) 

I  (2.57) 

Equation  (2.57)  is  the  dispersion  relation  for  the  3en-Wyller 
generalization.  Itr  simplification  and  solution  for  leads 
to  (see  Appendix  ft) 

i  A  +  B* in2<P  t  J  Bxi<nv(p- 

(A'lf)  -  -  J>+  ~  (2.58) 

where 

A  =  2(z  (eT  f  fsr)  (2.59) 

B  -  £fn  +•  6n? )  +  6jt 
c  r  ^  *x 

D  =  ^1 

£  r  i*£D  (2.63) 


(2.60) 

(2.61) 

(2.62) 
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Equation  (2.58)  is  the  Sen-Vyller  magneto-ionic  formula  and  will  be 
so  referred  to  in  what  follows. 


F  -  Ionospheric 


Licatior 


In  the  ionosphere  only  certain  kinds  of  velocity  dependence 
of  the  collision  frequency  are  of  importance.  For  example,  it  was 
indicated  at  the  beginning  of  the  chapter  that  Phelps  and  Pack  had 
determined  that  in  nitrogen 


v-  v  rr'x'* *  =  V  L 

y'yr»  1K.T  ™ 


(2.64) 


Since  nitrogen  is  the  most  abundant  gas  in  the  ionosphere  at  least 

at  the  lower  heights,  it  is  useful  to  compute  6Z  6^  €q j  for 

*  -r 

the  ionosphere.  If  <<■  kJ  ,  i.e.  if  the  field  applied  is  much 
less  than  the  thermal  energy,  then  by  equation  (2.16),  (2.18) 


=  a 


(2.6 5) 


f  (O) 

i.e.  x.  is  Maxwellian.  Therefore, 


,(»)  l3/x 

=  (anKT  /  e 


(2.66) 


The  handling  of  the  integrals  and  -T*  may  be  made  more 

convenient  if  they  are  written  in  their  complex  form,  keeping  ir>. 
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mind,  however,  that  it  is  the  real  part  the  one  that  is  of  interest. 
Their  complex  forms  are: 


j,  =  ciai  ((</,+  fc)<r*Jr* 

(2.6?) 

j,.- 

(2.68) 

J,=  e^f^+rK^ 

(2.69) 

It  is  easy  to  see  that  their  real  parts  are  in  accordance  with 
(2.36),  (2,37),  and  (2.38).  The  procedure  for  the  calculation  of 
(,  tir  and  will  be  illustrated  only  for  es  ,  the  rest 

I  j  U  — 

are  obtained  in  the  same  way* 

If  o'i  and  /3X  are  substituted  in  equation  (2.67)  then 


lut  rfzjL  v  '7P  '  —  I 

^  ,  -*  -  /a&k1, 

h  =  ,  51?*  "  V  KT  /  A 


(2.70) 


If  this 


but  since 

substitution  is  made  in  equation  (2.70)  and  some  simplification  is 
done,  then 


M 


S/2. 

L<*x^ 


,Vx  1 

LU3i  * _ 

+  J  (ir/A 


(2,71) 
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Therefore 


_  ,  Jiitc  trf  ,  ,  *  /■  ~  —  ]c_L  A 

*t:  ‘  ~£S  *  /4VV^  y*ft-*u*J~i7** i 


(2,72) 


If  some  rearrangements  axe  Bade*  then 


,  e&mh***  a  fr “^Cr, a  , 79, 


From  Burke  and  Hara  (Ref  23)  it  nay  be  seen  that  a  certain  kind  of 
integral*  called  the  semi-conductor  integrals  ara  defined  very  similarly 
to  those  in  equations  <2«?3)t  In  particular  the  (■»  semi-conductor 
integrals  are  defined  thus* 

/*?  ,  P  -(■ 

— —  /e 

r  f>!-0  fVr* 


(2.74) 


Therefore » 


✓ 

*-r  ' 


r  /- 


l'-.  *?4  (£)-&&  &  J  (2.75) 

ii  ^  2(*»> 


J/-* 

/ft 


or 


/-  £  As  A*  )  C2.7« 


-  /-  A  -  tb 


(2.77) 
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~7  ~~~ 


Where  &  -  —7-  (>u)  /  _  Csy  (£  )  *  Similarly 

elL  -  £U-J)+  f  (c-e)  (2.78) 

=£*~i(cH *.)]+  c[l>-±(-f+J)]  (2.79) 


Where 

c 


u)/(eO'-i) 


(2.80) 


*  •* 
■T^A, 


(2.8i> 


£  = 


(2,82) 


s(o£ 


(2.83) 


These  constants  a,  b,  c,  d,  e»  f  are  calculated  once  the  profiles 
of  collision  frequency  and  electron  density  are  known  and  therefore 
dexermine  (jJ  tjj >  ,  which  in  turn  will  determine  A,  B,  C,  D.  and  S 

in  the  Sen-Wyller  Magneto-ionic  formula.  So  then  M  f  ^  can  be 
determined  at  any  point  in  the  ionosphere. 
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A  Note  on  Sfakarof sky’s  Generalization 

la  this  Chapter,  the  classical  Apple ton-Har tree  theory  was  modified 
to  include  a  collision  frequency  which  was  dependent  on  velocity* 
Specifically,  the  results  were  specialized  to  the  ionosphere  by  using 
experimental  results  by  Phelps  and  Pack  which  shoved  that  for  nitrogen, 
the  collision  frequency  V  was  proportional  to  the  square  of  the  speed 
the  electrons.  The  assumption  of  a  slightly  ionised  gas  was  retained 
throughout  the  Sen-Wyller  generalization,  but  since  the  main  interest 
that  this  theory  has  for  the  author  is  its  application  to  the  ionosphere, 
it  becomes  clear  that  for  the  F-layer  and  above  where  the  number  of  ions 
has  increased  considerably,  a  slightly  ionized  gas  is  not  an  adequate 
model  of  the  ionosphere  for  these  heights.  Hence,  am  approach  which  could 
consider  any  degree  of  ionization  in  the  ionosphere  would  make  the  ideal 
approach.  This  approach  was  undertaken  by  Sbkarofsky.  Here  the  author 
outlines  Shkarof sky  *  3  approach,  shows  that  certain  approximations  ^ay 
be  used  and  justifies  a  fined  approach  for  the  ionospheric  absorption 
calculations. 

As  presented  in  the  first  two  chapters  any  generalization  of  the 
Appleton-Hartree  equation  must  undertake  to  generalize  the  conductivity 
tensor.  This  is  therefore  the  first  step  in  Shkarofsky's  generalization, 
the  generalized  conductivity  tensor  obtained  by  him  is  (Ref  5^)* 


£  r 


1+  c 

t  (b  #  -  c  ) 


c(6-c) 
b  +  C 
O 


a 


QgiggggSgEMBWXJli-fM  » «»**»**  •»  . .  ’  **  - 
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Where  cl  ~  J^Z.1  /  _ _ 

to*  +  <*>*o 


2c  ^ 


/Ye3  I 


m  /  <fo/s)A+ 


zb*  ^ 


Nez  / 


h*  <$>J-  +t(Os)t>_ 

where  is  the  total  average  collision  frequency  of  the  electron 

and,  y  -  and  ^  v  _  are  correction  factors  to  account  for 
the  variation  with  velocity  of  the  electron  -  neutral  particle 
collision  frequency  and  for  electron-ion  and  electron-electron  effects 
From  the  relationship 

(t  -  r^r  *■  i 

between  the  dielectric  tensor  and  the  conductivity  tensor  it  is  easy 
to  obtain  the  former  and  if  it  is  substituted  in  the  dispersion 
relation  it  is  a  cutter  of  algebraic  manipulation  to  obtain  a  new 
equation  for  the  refractive  index*  The  equation,  written  in  full,  is: 


/'iX(ze,  +i$.)  7 


-  izj  +{(-  +Z%  *  , 

*C*.-X-y.z)  ' 

K - ttk-  T7pj - J } 
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where  ^^4?  are  the  same  parameters  used  in  the  first  chapter;  and 
the  8  and  *</  functions  are  defined  as  follows:. 

ze,  =2j0~j-+j. 

29*  r  J*  -j- 

2S3  =  -  Y(Af  -A.J 

20,  -J.U.tiJ  t/i,(j.tfr)-ijJ.  -*jA.*n-jA*jJ>.-ij.ir  tij*!,.) 

--  fr  fr  tj.fr  -  2J.fr 

z4,  t>.  *L  +  Hi.-  A-) 

2  4*.-  fr  tfr 

*r%  *  4  (*.  +  ij-  >1.^+  u.(K  -  k)  td.fr  r 

2/4,  -  +  Y(fr  t 


and  f.j  h,  are  the  values  of  j}  /  for  argument  j  J+j  for 

argument  ,  and  J~}  f>_  for  argument  . 

The  behavior  of  these  functions  should  be  explained*  For  a  given 
value  of  <****^/<y £«>  (the  ratio  of  collisions  with  atoms  over 
collisions  with  ions,  which  indicates  the  degree  of  ionization  of  the 
gas;  for  high  values  the  gas  is  slightly  ionized,  low  values  tend  to 
indicate  fully  ionized  gases)  the  arguments  are  made  to  vary;  j 
and  h  are  then  seen  to  be  roonotonic;  for  small  values  of  , 

y  increases  from  a  limiting  value  less  than  unity  to  unity  at  large 


I 

i 


i 

4 
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values  of  ^U~SV<^'>  t  while  >»  decreases  from  a  limiting  value 
greater  than  unity  to  unity. 

The  behavior  with  respect  to  is  somewhat  more 

complicated.  It  is  related  to  the  velocity  dependence  of  the  electron- 
neutral  particle  collision  frequency 

H,  -  cir^ 

where  o'  is  the  velocity  of  the  electrons;  if  r* o  ,  J  and  k  have 
a  monotonic  variation  with  ^ if  j  and  h  ,  are 

unity  for  large  values  of  4 ^ >  ;  and  if  h>(  ,  j  and  k  vary 
from  the  values  for  the  slightly  ionized  gas  to  approximately  unity, 
at  between  1  and  10,  and  then  j  and  k  deviate  more  from 

unity  as  decreases,  until  they  approach  the  limits  of  the 

completely  ionized  gas. 

Hence,  if  4J^s4»dpand  to>S  ,  then  all  the  j  and  k  functions 
approach  1  and  the  usual  Appleton-Hartree  equations  result  with  ~ 
how  the  conditions  cited  here  can  apply  to  the  F  region  of  the  ionosphere. 
Shlcarofs&y  has  determined  as  a  good  estimate  that  if 

S  > 

then  investigation  of  the  6  and  v£  curves  reveal  that  it  is  safe  to  use 
the  ordinary  Appleton-Hartree  equation  provided  that 

if  .5  °-t  too  i  <  z. 


6ij 


f»-ts  Antt  //  h  « 

>jr/rzi/  ojj-A 


In  the  present  paper  all  the  frequencies  used  are  h.gher  than  2*5 
Mc/s  so  that  the  first  condition  is  always  satisfied.  Hence,  we 
can  combine  the  first  two  approaches  to  describe  a  variable  state 
of  ionization  in  the  ionosphere:  below  the  heights  at  which  the 
aforementioned  conditions  begin  to  apply  the  Sen-iv'yller  theory 
should  provide  a  valid  description  of  the  ionosphere  and  at  heights 
where  the  condition  applies,  the  ordinary  Apple toa-Hartree  equation 
might  be  used. 
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IV.  NUMERICAL  APPLICATIONS 


In  this  Chapter  the  procedure  followed  will  he  described.  For 
the  sake  of  illustration,  the  example  of  analytical  profiles  for 
the  electron  density  /VU)  and  collision  frequency  ?(£■)  will  be 
used  to  compute  absorption  and  the  frequency  dependence  of  it. 

The  particular  profile  of  /] f(k)  used  as  an  example  is  the  Chapman- 
layer  profile  although  there  exists  some  other  profiles  which  lend 
themselves  to  easy  computations.  The  Chapman-layer  is  however,  the 
more  realistic  and  useful.  For  a  discussion  of  the  other  profiles, 
Sudden  (Ref  12:  Ch  10)  may  be  consulted. 

After  the  examination  of  this  particular  profile,  the  profiles 
used  in  the  present  calculations  will  be  introduced  and  some  words  as 
to  the  manner  of  their  selection  will  be  set  down,  A  discussion  of 
several  difficulties  which  are  encountered  in  any  calculation  of  this 
type  will  be  given.  Also,  possible  errors  and  systematic,  sources  of 
inaccuracies  Jill  be  discussed  along  with  certain  effects  which  are 
neglected  in  the  computations. 

A  -  The  General  Line 

As  mentioned  at  the  end  of  the  first  chapter,  the  fundamental 


equation  in  these  calculations  is 


MU)  - 


(3.1) 


Where  A  is  the  absorption  of  the  wave  in  db,  ds  the  path  of  the  wave 
and  /r  the  absorption  coefficient.  For  each  sign  in  the  Appleton- 
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Hartree  or  the  Sen-Wyller  equations  a  different  ^  can  be  determined. 
These  will  be  designated  by  and  ,  the  first  for  the 

ordinary  ray  and  the  second  for  the  extr odinary .  The  differential 
absorption  is  defined  as: 

$  -  A*  "At  (3.2) 

Therefore 

t>  -  4-33  -{.JSs  (3.3) 

One  of  the  objectives  of  this  paper  is  to  determine  the  frequency 

dependence  of  A* ^  A ,  and  D  •  Two  things  are  necessary  for 

the  calculation  of  these  integrals::  know  the  path  of  the  ny  and  faa 

as  a  function  of  the  path.  Since  the  usual  experiments  restrict 

themselves  to  the  vertical  direction,  the  path  and  the  height  may 

% 

be  said  to  be  equivalent,  This  point  needs  a  little  more  discussion, 
for  it  is  not  q<.  j.te  so  simple  as  that,  but  it  will  be  posponed  until 
later  after  the  general  direction  of  the  calculations  has  been  outlined. 
The  other,  the  determination  of  fa,  as  an  analytic  function  of  height, 
is  something  that  may  be  done  only  in  the  simplest  of  cases.  This 
will  be  illustrated  in  the  example  included  in  tuis  chapter.  It  is 
convenient  to  express  ^  as  a  function  of  the  real  and  imaginary  parts 
of  the  complex  index  of  refraction  which  was  the  quantity  derived  in 
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in  the  first  two  chapters,  Ibis  has  already  being  done.  The  equations 
ares 


*  =  7z  f**' Jr* +7* 


(1.113) 


Jr2  7t 


(1.114) 


where  R  is  the  real  part  of  W  .  and  I  the  imaginary.  So  that 

ft  -  4  33 J ^  (X,  l)jl  t3‘4) 

A  and  ■  are,  of  course,  functions  of  height. 

Another  definition  is  of  usefulness  at  this  point.  In  many 
instances  A  may  be  written  as 

f.l  ifKNMJJ.  (3-5) 

Where  JC  -  ir(v,w)  independent  since  it  represents  the  absorption  per 
electron  per  kilometer  of  path  as  a  function  of  height.  When  this 
separation  is  possible  considerable  labor  is  saved  since  K  can  be 
calculated  independently  of  the  electron  density  profile. 

Lastly,  the  limits  of  the  integrals  are  taken  according  to  the 
region  of  interest.  For  the  total  traversal  of  the  ionosphere  the 
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the  upper  limit  ia  taken  according  to  the  estimates  of  the  remaining 
ionization  higher  up*  A  cutoff  point  is  selected  on  the  basis  of 
the  estimated  error  which  would  result  if  the  upper  limit  were 
chosen  at  a  certain  particular  height* 

B.  -  Chaaaan-Layar  absorption  for  the  31.  Approximation  in  the  A-H  Theory 

In  this  exarjple  (fief  27)  some  features  of  the  absorption  calculations 
are  brought  out*  It  serves  as  a  guide  for  the  more  complicated  calculations 
which  follow* 

For  the  3JL  approximation  in  the  Appleton-Hartree  equation* 

2  _X_ _ 

"  i-iZ  £  K.  0.6) 

Where  the  +  sign  denotes  the  ordinary  ray*  the  -  sign  the  extraordinary 
ray.  Equation  (3*6)  may  be  put  in  the  form 

M*  -  l-  CM  r 

xOsfJ  t  Z±- 

If  for  the  particular  layer  being  examined*  then  ^  -v  / 

(the  absorption  undergone  when  m*1  is  called  non-deviative  absorption) 
and  \  •  Therefore, 

ZX  (3.8) 

f r  (i±K.)y+zi 
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In  taxes  of  frequencies  the  absorption  coefficient  is* 

f 1  <[o>tsr*  p‘]  (3.9) 

a  tie i 

If  it  is  recalled  that  ^  '  $  fa  and  if  tots  -  6)t  (the  effective 
operating  frequency),  then 


The  specific  absorption  is 


e  *_  JL _ _ 


(3.10) 


(3.H) 


which  is  a  function  of  height  because  and  S  functions  of  height. 
It  should  be  stressed  that  in  this  approximation  the  collision 
frequency  is  not  negligible  when  compared  to  the  effective  operating 
frequency.  Because  of  this  reason  this  equation  could  be  applied  to 
D-region  absorption  when,  the  other  conditions  of  high  frequency 
(compared  to  t'ne  critical  frequency)  obtain.  It  may  also  be  seen 
in  what  way  profiles  of  electron  density  and  collision  frequency  are 
needed  for  the  computations  of  absorption.  For  a  realistic 

* 

atmosphere  these  profiles  have  no  analytical  expression,  a3  mentioned 
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in  the  introduction. 

Some  useful  computations  may  be  performed  if  a  particular  ionized 
layer  is  assumed  to  be  of  the  Chapman -Layer  type#  i.e. 


Nfc)*  A£ex/>  t('“m 


(3*12) 


Where  'K  is  the  peak  electron  density  when  the  sun  is  at  the  zenith* 

A- A, 

^  is  the  zenith  angle  and  Z  ~  —jf  .  The  collision  frequency 

.  -z 

is  assumed  to  hare  the  profile  r»  C  where  K  is  the  collision 

frequency  at  A  -  m 

If  /V(A)  and  y(A)  are  substituted  in  equation  (3*10)  then 

e  *  K  ~  f  '-j2  -  sec.fC  ~Z)&  £ 


T“ 


4,o»  y/e'2*] 


(3.13) 


Since  y  is  now  explicitly  a  function  of  *  there  should  be  a 


change  of  variable  free  k  to  Z  in  equation  (3*4)*  Since  z~  ~jj 
then  c/A :  Afdz  .  Therefore 


(3.14) 


A- A, 


,  /  jt  .  f  efy>(Z-ljz  -Zjecgre~£)  , 

Ai  =  4.33 j -  ,.33  £oCM  J  fQx+  y/’e'** J 
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The  operating  frequency  u>c  is  a  constant  if  the  gyrofrequency  is 
constant*  In  the  following  steps  it  will  be  assumed  constant* 

If  the  substitution  y  •=  / j^Sec  tj  i't  la  made  the  integral 

J  1 

can  be  simplified  considerably*  A  is  then 


4.55  (  . 


C  >+} 


Ac  I  f  ~~  VTcxpO)  cos 


3s. 

1 LI 


(y  v/ * ) 


*y 


where 


/3  - 


£JeSecy 

*K 


or 


-  +.33 


c  cP. 


:( 

J  cos 


■>{(/2'Vo)  Jy 


<3.1 5) 


Since  Z  goes  from  -»o  to  &o  %  y  goes  from  o*  to  o  • 
Therefore  equation  (3*15)  can  be  written  as 


A\  r  4.33 


(3-16) 


Where 


If 


(3.17) 
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s  */>'/¥  kW 


(3.18) 


Where  Cf>&)  is  defined  as  in  the  second  chapter*  It  is  one  of 
the  semi-conductor  integrals.  These  have  been  tabulated  by  Burke 
and  Kara-  (Kef  23).  For  some  of  the  £  integrals  (the  ones  with 
p  =  1.5  and  2,5)  polynominal  approximations  have  been  derived 
empirically  by  means  of  computers.  The  author,  however,  has  not 
seen  one  for  p  =  *5*  The  usual  thing  in  this  case  is  to  express 
(p)  as  a  function  of  the  Fresnel  integrals  (for  details  see 
Jaeger  (Ref  2?)  and  to  evaluate  these  numerically.  The  advantage  of 
the  polynomial  expressions,  apart  from  their  ease  in  computation,  is 
that  they  allow  the  frequency  dependence  of  absorption  to  be  examined 
more  clearly.  At  first  sight 


A*  Zl 


O.19) 


Where  is  the  effective  operating  frequency,  but  in  reality,  of 


course 1 


/!  <*'  — ”  JO  (&t ) 

(t)x 


(3- 20) 


According  to  Jaeger* s -tabulation  of  ■?(/&) 


i(F>  - 


constant 


(3.21) 
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when 


.  Therefore*  when 


a>t  > 


wA0 

J*c{ 


t  then 


J  ,Ti  (3.19) 

ut 

but  before  that  condition  is  met  equation  (3*20)  is  the  one  that 

is  valid*  Interestingly  enough,  in  terms  of  ^  (p)  ,  equation  (3*20) 

is 


A-r  (*>S  C'/x  fat)) 


(3.22) 


A* <  C'/J*0*)  (3,23) 

The  differential  absorption  would  be  proportional  to  the  difference 
of  the  integrals  for  the  corresponding  • 

For  other  profiles,  it  is  convenient  to  calculate  MCA)  and 
then  multiply  it  by  the  particular  MU)  at  the  height  involved* 

From  (3*11) 


(3.24) 


Differentiation  shows  that  this  function  K  has  a  maximum  when 
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Sinew  u)^  ~  cofs  »  equation  (3»25)  reveals  that  at  the  heights  at  vhich 
K  is  maximum  for  the  two  magnetoionic  components* 


Uf,  =  y,-s 


(3.2 6) 


r  K  +S 


(3* 27) 


which  means  that  H  S  y2  •  This  implies  that  the  maximum  of  k  is 
higher  for  the  extraordinary  ray  then  for  the  ordinary  since  it  requires 
a  lower  collition  frequency  than  the  ordinary  ray  to  satisfy  condition 
(3.26) 

Another  very  interesting  fact  can  be  obtained  from  equation  (3. 10). 
If  the  absorption  coefficients  for  the  extraordinary  ordinary  rays 
are  denoted  by  and  p  respectively  by  denoting  the  respective 
operating  frequencies  by  cOc*  and  cJfm  ,  then 

4L  .  KUtl 

f.  o*28) 

which  is  independent  of  •  When  /^K  and  po  are  plotted  on  logarithmic 
scale 


(3.29) 


and,  furthermore,  when  then  /oj  f*  ib  »®srly  constant, 
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It  should  be  stressed  that  the  previous  equations  are  strictly  valid 
only  under  the  specific  conditions  under  which  they  were  derived.  In 
the  ionosphere  they  will  not  always  be  valid  but  there  are  certain 
ranges  of  heights  in  which  they  can  be  used  to  simplify  the  analysis. 


Special  Cases  in  the  Sen-Wy] 


For  longitudinal  propagation,  the  Sen-'.v'yller  magneto-ionic  formula 


reduces  to 


»>.  //- 


(3.30) 


If  the  frequencies  O  are  high,  then  /a./  ,  and 


„  s<J,1  p  /  J 

AT  -  *  (»**  ) 

f-  zz f*  (^J 


(3.31) 


(3.32) 


/js iJL£  ) 

The  specific  absorption  function  <  ~  fjZjz  /  *  wJlich 

is  again  independent  of  N.  Another  fact  which  should  be  stressed  is 


that 


jfi  _  ) 

!-  e*  (*gj 


(3.33) 


76 


sp/ph/65-i 


or 


in  terns  of  the  effective  operating  frequencies* 


So  that  Judging  from  Burke  and  Hara  (fief  23)  for  high  altitudes  where 
is  very  low,  4r  should  be  nearly  a  constant.  (See  Fig  2).  There 

r 

are  other  valuable  approximations  which  can  be  made  but  they  are 
not  useful  in  this  case*  For  a  list  of  this  approximations  Hatcliffe 
(Ref  45)  can  be  consulted  along  with  the  original  papers* 


D  -  Selection  of  Profiles 


Perhaps  the  previous  sections  have  given  an  indication  of  two  ways 
that  one  can  proceed  to  compute  absorption  in  the  ionosphere*  One  can 
either  select  several  profiles  for  electron  density  and  a  collision 

i 

frequency  profile  or  one  can  compute  the  specific  absorption  function 

,  I 

and  then  multiply  by  the  particular  profile*  The  second  procedure  has  | 

the  advantage  that  the  once  the  function  has  been  obtained  one  can  | 

multiply  it  by  any  desired  profile  and  obtain  a  quick  idea  of  the  1 

variation  of  the  absorption  coefficient  but  then,  for  complicated  f 

ionospheric  conditions  there  is  no  guarantee  that  K  is  independent  of  •  j 

The  previous  sections  indicate  that  K  is  certainly  independent  of  N  for  j 

I 

longitudinal  propagation  but  one  has  to  be  careful  that  this  condition  j 

I 

obtains  to  use  it  for  any  particular  interpretation*  | 

5 

r 

Both  procedures  were  followed  in  this  paper.  First  the  selection  J 

of  the  profiles  used  will  be  described,  and  in  a  later  section  some  | 

1 

mention  will  be  made  of  the  other  procedure*  | 

The  profiles  of  electron  density  and  collision  frequency  selected  1 


3- 
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for  computation  are  necessarily  of  a  representative  nature*  Among 
the  reasons  for  this  are  the  inaccuracies  present  in  the  actual 
measurements  of  these  parameters  in  the  ionosphere*  the  complex  set 
of  conditions  that  determines  these  parameters*  and  daily*  monthly, 
and  other  variations  which  result  in  a  large  spread  of  values  for  ^ 
and  y •  Furthermore,  in  addition  to  being  merely  representative  they 
are  also  necessarily  extrapolated  in  certain  cases;  the  author*  for 
example*  knows  of  no  single  experiment  which  has  measured  ^  through¬ 
out  the  whole  ionosphere.  Therefore  several  groups  of  data  from 
different  experimenters  were  tied  together  by  appropriate  extrapolation 
to  form  a  composite  profile  which  is  assumed  representative  throughout 
the  ionosphere. 

For  the  daytime  profile*  the  lower  heights  were  covered  by  a 
profile  derived  by  Thrane  and  Piggott*  at  Kjeller,  Norway  on  the 
basis  of  cross-modulation  experiments.  The  profile  expeci&lly  applies 
to  noon-eq.uinax  conditions  at  sunspot  maximum.  Values  from  this  profile 
were  plotted  up  to  a  height  of  100  km.  From  heights  between  100  and 
200  a.  profile  by  A.  K.  Paul  and  J.  ¥.  Wright  (Ref  43)  obtained  at 
Eglin  Air  Force  Base*  Florida,  was  used.  These  authors  have  calculated 
several  profiles  for  different  hours  of  the  day.  The  one  st»lectel  was 
a  noon  profile  which  most  smoothly  fitted  the  one  at  Kjeller;  the  range 

of  electron  density  in  the  heights  where  there  was  overlapping  was 

JfJ 

smoothed  out  by  a.  very  small  change  in  at  the  end  point  in  the 

*  Unpublished  puper* 


““f  *■  V  yr  _  >  ”Sg  f 4  ~~ 


k 


| 
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curve  of  Thraae.  The  rest  of  the  daytime  profile  is  obtained  free 
a  profile  by  Shkarofsky  (Ref  5°)  up  to  a  height  of  34^ V.m  •  Beyond 
340  a> linear  decrease  of  N  is  assumed.  This  assumption  is 
corroborated  by  topside  sounder  data  for  dip  angles  of  approximately 
80°  (See  LeGalley)  (Ref  33*192). 

The  results  of  two  authors  were  used  for  collision  frequency 
during  the  day:  those  of  Nicholet  (Ref  39)  and  Shkarofsky  (Ref  5°)* 

Those  of  Nicolet  were  used  for  the  lower  heists  (up  to  about  100  k»*) 
and  for  the  rest  of  the  heights  those  values  given  by  Shkarofsky. 

A  word  of  caution  should  be  added:,  the  curves  given  for  collision 
frequency  of  the  electrons  include  collisions  with  ions  which  became 
of  importance  in  the  F-region;  the  Appleton-Hartree  and  the  Sen-Wyller 
formula,  however,  include  only  electron-neutral  collisions;  since  in 
the  lower  ionosphere  the  electron-ion  collision  frequency  is  negligible 
compared  to  the  electron-neutral  collision  frequency,  the  Appletoa- 
Hartree  and  the  Sen-Wyller  equations  are  valid  (with  reference  to 
collision  frequency  considerations)  for  the  lower  ionosphere.  As 
the  neight  increases,  they  become  less  valid.  Obviously,  the  variations 
of  %  and  >4*,  with  height  are  partly  a  result  of  the  varying  degreo 
of  ionization  of  the  ionosphere  as  the  height  varies*  A  derivation 
of  the  magneto-ionic  formulas  which  would  recognize  this  fact  would  be 
preferable.  This  approach  has  been  tried  by  Shkarofsky  (Ref  3°)  (See 
note  included  in  this  paper). 

The  night-time  N(k)  profile  was  selected  from  an  KBS  Technical 
Note  (Ref  56).  The  profile  was  derived  for  the  month  of  January  i960 
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at  Puerto  Rico.  The  values  of  Hlk)  are  provided  in  the  Note  in 
tabular  form  for  the  24  hours  of  the  day.  Of  these,  the  0400  o’clock 
profile  was  selected  as  it  was  that  one  with  the  smallest  peak  in 
electron  density  and  so  deviated  the  most  from  a.  daytime  profile.  The 
tabular  values  were  plotted  and  joined  by  straight  line  segments.  V/ith 
this  approximation  in  mind,  a  new  table  of  values  was  made  for  h  a> 
which  filled  the  gaps  the  other  had* 

The  night-time  collision  frequency  was  mere  difficult  to  get. 
Night-time  collision  frequency  profiles  are,  in  fact,  not  readily 
available  in  the  literature,  a  profile  had  to  be  assumed  which  would 
provide  a  reasonable  basis  for  night-time  calculations  of  the  index  of 
refraction,  For  heights  below  200  km  the  same  exponential  decrease 
was  used,  above  200  Km  the  night-time  /$) profile  vras  used  to  recalculate 
collision  frequency  values  according  to  formulas  given  by  Shkarotsky 
(Ref  5°)*  The  two  curves  were  joined  by  extrapolation  following,  to 
a  certain  extent,  the  shape  of  the  day  time  profile# 

as  can  be  seen,  the  selection  of  these  models  was  to  a  certain 
extent  a.  oitrary.  They  serve,  however,  a  purpose,  i.e.  to  bring  out  of 
the  formulas  certain  regularities  or  features  which  are  Independent  of 
the  model  chosen.  The  conclusions  are  then  relative  things,  coaparisions, 
frequency  dependence  relations,  etc* 

E  -  The  Calculations 

The  calculations  are  based  on  the  following  two  formulas  derived 
in  the  first  and  second  chapters,  the  Appleton-Fartree  formula: 
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X 


t  [XL — r„  /  V 

t /*{/-/ -tZ)* 


and  the  Ssn-Wyller  mgneto-ionic  formula: 

2_  Bsin  V  X  -cjcm*0 

P  -f-  iz  s/y>*  (f 

It  was  previously  stated  that  if  M 1  was  written  as  X-el  i/here 
ft  and  -/  are  the  real  and  imaginary  part  a  of  //^then 

m*7iJk  +  /r*  +  iz 

?  "  J ~  ft  4-  Jit*  +- 1 z 

ft  and  X  *«■*  ia  general  complicated  functions  of  /V  and  y  through 
which  they  depend  on  the  height,  except,  of  course*  in  cases  as  simple 
as  the  one  presented  in  the  beginning  sections  of  this  chapter.  Numerical 
calculations  had,  therefore,  to  be  undertaken.  Two  different  programs 

were  used  (one  for  the  Appleton-Hartrea  and  the  other  for  the  Sen-Wyller 

.v 

formula)  to  separate  M  into  its  real  and  imaginary  part.  Given  a 
certain  geographical  location  the  angle  of  propagation  between  a  vertically 
incident  wave  and  the  magnetic  field  can  be  determined,  which,  if  there 
is  no  deviation,  is  constant  for  a  dipole  field  (See  appendix  D);  K-  ” 
can  be  determined  at  any  height  if  the  variation  of  gyrofrequency  for 

a  dipole  field  is  included,  i*e« 

r  \3 

/  _  I 

^  +.  k  i 
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Whero  s,  is  the  electron  gyrofrequency  at  ground  level,  ro  the 

h  j 

radius  of  the  earth  and  for  the  height.^  ,  X  and  /-  can  be  determined, 
f^om  the  collision  frequency  and  electron  density  profiles  respectively , 
as  functions  of  height-*  For  each  height  then,  XL  and  7  c3*1  he  computed. 

For  the  Sen-Vr'ylier  formula  the  ease  is  true,  although  previous 
tabulation  of  the  semi-conductor  integrals  (°  (%)  .  JL  j  |b-iL  c/(-  is 
necessary.  These  can  be  tabulated  from  formulas  given  by  Burke  and 
Kara  (r*sf  tf3\5~6)* 

With  the  purpose  of  allowing  the  frequency  dependence  of  certain 
quantities  to  be  brought  out  more  explicitly,  a  vide  range  of  frequencies 
was  used  for  computation.  These  frequencies  vers  2,  5*  3*  5»  ?*  10 »  14» 
20*  itO,  80  .  Some  of  these  were  examined  under  both  day-time  and 

night-time  conditions. 

Once  ** ,  v  are  computed  as  functions  of  height  then  numerical 
integrations  ssy  be  performed  os  them  tc  calculate  quantities  defined 
at  the  beginning  of  this  chapter. 

Finally,  the  geographical  location  used  (Sagamore  Hill,  Mass)  is 
consistent  with  Q  ~  73°  and  Sp  =  1.57  ,  and  YT  V*.  are  computed 

from  these  for  any  frequency  and  height, 

F  -  Wot- a  on  Possible  Sources  of  Inaccuracies 

The  calculations  done  by  the  simple  method  explained  above,  involve 

no  approximations  so  far  as  ^(k)  and  /-* ( 7j  are  concerned.  But  it 

i 

will  be  recalled  that  to  evaluate  the  integrals  .the  approximation  was 
«Qi*de  tnat  the  path  of  the  wave  could  be  set  equal  to  the  height,  this 
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approximation  is  used  throughout  the  literature  and  is  valid  for 
high-frequency  waves.  Without  going  into  much  detail,  however,  it 
may  be  mentioned  that  there  occurs  always  some  deviation  in  the  path 
of  the  ray.  This  deviation  becomes  particularly  important  when  m -v  ' 
becomes  much  less  than  unity  (i.e.  for  low  frequencies).  Since  the 
interest  of  the  present  paper  is  on  waves  which  penetrate  the 
ionosphere  (see  lut.),  it  should  be  stressed  that  the  carves  computed 
for  2,5  ** %  s  5  yv'//5  »  7  ^ V$  can  not  be  integrated 

throughout  the  whole  ionosphere  since  these  are  reflected  somewhere 
in  the  ionosphere.  It  was  also  assumed  in  these  calculations  that 
the  field  of  the  earth  was  a  dipole  field.  This  is  only  approximately 
true,  a  first  order  approximation  as  extensive  investigations  (itef  33: 

Ch.  9)  of  the  geomagnetic  field  show.  A  third  factor  which  hinders 
the  accuracy  of  the  computations  is  the  accuracy  of  the  electron 
density  and  collision  frequency  profiles,  which  are  determined  by  several 
different  methods  under  varying  conditions.  The  best  one  can  hope  to  do 
is  to  use  profiles  for  average  representative  conditions  such  as  mean  noon 
electron  density  profile,  etc. 

.x  fourth  factor  which  is  implicit  in  the  above  calculations  and 
which  also  effects  the  accuracy  is  the  Sea-syller  approach  itself.  As 
was  already  mentioned,  the  Sea»-«yller  approach  is  strictly  valid  for  a 
slightly  ionised  Lorentz  gas.  The  ionosphere  will  satisfy  this  condition 
at  the  lower  heights  but  probably  not  at  the  F-regioa.  3ifcarofsky*s 
generalization  considers  all  degrees  of  ionization.  Unfortunately  certain 
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V.  RESULTING  HEIGHT  PROFILES  OF  U  t  {  t&D  ASSOCIATED  RESULT 

The  results  of  the  different  computations  performed  are  presented 
in  this  Chapter.  They  represent,  in  the  first  place,  height  profiles 
for  the  integrands  for  ordinary  and  extraordinary  absorption.  Also 
some  profiles  for  My  have  been  presented.  The  computed  profiles 
were  graphed  automatically,  and  are  presented  in  Appendix  F,  A  stands 
for  N(L)  t  8  for  y(^)  t  U+  for  »  snd.  o  ;  the  log  to 

the  base  10  of  each  of  them  has  been  plotted  against  height.  The  main 
purpose  for  including  them  is  because  they  are  necessary  for  the  calculation 
of  the  absorption  integrals,  but  the  h'  ght  profiles  of  ^  ,  and  %  are 
of  interest  in  themselves.  The  interpretation  of  the  shape  of  these 
curves  is  not  a  simple  matter,  and  it,  being  not  strictly  germane  to 
the  present  calculations,  has  not  been  done  extensively,  but  in  the 
main  has  been  restricted  to  several  general  statements. 

In  the  second  place,  the  curves  have  been  also  grouped  into 
significant  pattern* to  perait  easier  discussion.  These  are  presented 
in  Appendix  G. 

Finally,  seme  curves  of  absorption  coefficients,  exploring  its 
proportionality  to  ,  if  it  exists,  have  been  presented  in  Appendix 

H. 


A  -  General  Description 


The  values  of  //  and  )  computed  by  using  the  Sen-Wyller 

program  were  plotted  as  functions  of  height  for  both  day  and  night-time 


conditions.  Not  all  of  the  possible  ^ K  have  been  plotted. 

The  frequencies  used  were  2,5.  5»  7.  10,  li»  and  20  .  See  Appendix  F, 
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Of  the  curves  plotted,  several  bread  features  may  be  noted.  For 
night-time  conditions,  and  for  frequencies  above  the  critical  frequency 
of  the  F-iayer,  (5.6  mc/s)  X+  and  X-  are  seen  to  follow,  in  general, 
the  snape  of  the  electron  density  profile.  It  is  also  seen  that  it  is 
the  extraordinary  ray  which  first  departs  from  this  simple  behaviour 
as  the  graphs  for  a+,  X-  at  5  mc/s  show:  while  X+  still  has  a  variation 
similar  to  the  one  pointed  out  above  while  X-  has  a  fairly  "thick"  and 
distinctive  peek.  For  i requeue ies  belo  the  critical  frequency  they 
snow  marned  behaviour  with  height,  the  behaviour  of  X-  being  the  more 
complicated. 

Correspondingly,  the  curves  forX+  show  a  distinct  behaviour  at 
the  frequency  changes.  For  frequencies  above  the  critical  frequency 
they  are  very  nearly  constant  (approximately  1)  but  as  the  frequency 
of  the  wave  decreases,  a  slight  depression  begins  to  appear  in  the 
curves  of  X+  near  jCO  km  until,  for  frequencies  as  low  as  2.5  mc/s 
this  depression  has  turned  itself  into  a  valley  o i  lcvij{+  values.  The 
curves  for  M  serve  to  indicate  regions  of  deviative  absorption. 

or.cier  ^ay-time  cordi  tions  the  fcenaviour  of  the  plotted  curves  is 
much  more  complicatec,  but  again  broad  features  may  be  noted.  For 
frequencies  aoove  tr.e  critical  frequency  of  the  ionosphere  as  determined 
from  the  day-time  profile,  (about  5  mc/s).  the  behaviour  of  X+  and  X- 
(with  the  exception  c-f  X-  for  1C  mc/s)  above  3^0  km  is  similar  to  that 
of  the  electron  density  profile.  Below  that  height  they  show  two 
smaller  peaks  at  low  heights.  For  frequencies  below  the  critical 
frequency,  the  behaviour  of  X^,  X-  is  much  too  complicated  to  be 
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described  by  words,  but  again  is  the  more  complicated  curve  and  is 
characterized  by  extremely  sharp  peaks  at  certain  heights. 

The  /{t  curves  show,  in  day-time  conditions:  they^>  for  the  hign 
frequencies  are  fairly  constant,  while  for  decreasing  frequencies  they 
tend  to  dip  at  the  region  of  maximum  electron  density, 

A  limited  comparison  of  these  features  may  be  found  in  Appendix  G, 
where  ^  ^ g  and  ytf,  have  been  plotted  together  for  frequencies 

near  the  critical  frequency  and  above  for  both  sets  of  data.  The 
reasons  for  using  only  these  is  that  the  integrations  for  absorption 
were  going  to  be  performed  only  for  frequencies  which  are  able  to 
penetrate  the  whole  ionosphere. 


B  -  Absorption  Per  Electron 

Another'  calculation  for  which  results  have  been  obtained  is  the 
investigation  of  the  relation  between  the  absorption  per  electron  and 
the  electron  density.  In  a  previous  chapter,  it  was  argued  that  if 
jr  KN  t  where  K  was  a  function  only  of  frequency  and  collision 
frequency,  then  the  calculation  of  absorption  became  all  the  easier.  For 
longitudinal  propagation,  this  was  shown  to  be  true  but  for  general 
directions  this  may  not  be  true.  Therefore,  a  start  was  undertaken 
towards  investigating  this  problem. 

In  particular  /°C J  were  calculated 

and  these  were  examined  to  see  whether  Q'(V,N)  was  equal  to  KN  t  i.e* 

l 

if  %  then  /)-  K  ,  30  that  for  any  fixed  frequency  /)N 

provided  KN  • 
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In  general*  one  may  formally  define  a  function  %  such  thatp  7t^  , 
where  ?[  ~  .  Under  certain  circumstances  should  be  identical 

to  <  ;  it  is  these  circumstances  which  were  started  to  be  investigated 
at  the  end  of  the  coop  period. 

Some  of  the  curves  obtained  can  be  seen  in  Appendix  H.  They  cover 
a  v-  Va  range  of  frequencies:  80,  40#  20,  10,  5»  3»  2.5  and  1*5  Yr'c/s 
tor  the  angle  of  propagation  0  =  1.2?  radians. 


C  -  Frequency  Dependence  of  Absorption 

It  has  been  indicated  previously  that  because  no  analytical  profiles 
existed  for  N  and  V  the  absorption  had  to  be  calculated  numerically. 

The  plots  presented  in  Appendix  F  are  for  ^  and  but  these  are  merely 
the  respective  integrands  for  ordinary  and  extraordinary  absorption 
divided  by  ~  •-  Since  the  present  paper  is  only  interested  in  the  absorp¬ 

tion  undergone  by  a  wav®  which  comes  from  outer  space,  the  absorption 
integrals  were  calculated  only  for  those  frequencies  which  could  penetrate 
the  ionosphere  under  the  assumed  conditions. 

These  integrations  were  done  using  what  amounts  essentially  to  the 
trapezoid  rule.  This  rule  of  integration,  according  to  books  in 
numerical  analysis,  gives,  for  sufficiently  small  size  of  the  interval 
K  of  integration,  that 


/'1/X A  fy  +  A  ft  * 


a 


(5»l) 


Where  £jo  -  tj  ) -  j ( K0]  .  In  this  paper  the  integration  was 

merely  approximately  by  a  sum;,  thus:; 


(5.2) 
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Where  the  sum  is  carried  between  the  appropriate  limits  and  d  h  -5  km. 

,  +  ft  . 

From  (5.1)  above,  it  may  be  seen  that  when  J*.  y “X  is 

extended  to  a  large  range  for  dx  then  Jjjfa  =  f  x  ' 

The  difference  between  this  and  (3*2)  is  merely  r  *£j3  which  in 

the  present  case  is  •£*  4*  £%  and  this  is  very  small  for  any  frequency 

2*  2. 

judging  by  the  profiles  presented  in  Appendices  F  and  G, 

The  trapezoidal  rule  itself  is  a  curde  approximation,  but  first, 
may  be  considered  to  be  smal‘1  enough,  this  is  strikingly 
demonstrated  by  the  graphs  in  Appendix  G,  and  second,  the  main  interest 
in  performing  these  integrations  is  in  bringing  cut  the  frequency 
dependence  of  absorption,  it  is  therefore  felt,  that  given  the  inaccuracies 
present  in  the  data  itself  any  more  specialized  rule  such  a3  Simpson's 
rule  will  not  bring  out  any  new  fact  about  the  frequency  dependence. 

It  was  estimated  that  the  contributions  to  the  integral  from  heights 
above  700  km  were  less  than  \%  of  the  total  integral,  hence  the  integrals 
were  carried  out  to  700  km. 

It  should,  at  this  point,  be  emphasized  that  these  integrals  include 
aeviative  absorption  that  might  be  present  since  no  particular  assumpt- 
tion  was  made  about  /<  ,  although  no  attempt  was  made  to  separate  the 

contributions  from  both  components,  and,  as  was  mentioned  previously, 
no  attempt  wac  made  either  to  estimate  rigorously  the  path  of  the  waves, 
something  which  constitutes  a  problem  unto  itself. 

As  a  matter  of  convenience  the  integral*  were  calculated  up  to  375 
km  from  375  km  upward  for  night- time  conditions.  Thi3  made  it 
possible  to  tali:  about  contributions  from  different  parts  of  the 
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ionosphere  to  the  total  absorption.  For  tne  day-time,  the  following 
segments  were  Integra  tec  separately:  up  to  1C5  km,  from  1C-5  to  200  km, 
2C0-310  km,  anc  310-70C  ran,  Fach  of  these  regions  includes  a  certain, 
ionospheric  region  of  interest  as  may  be  seen,  Irom  the  plots  in  Appendix 
G,  The  results  are  presented  below,  where  Aq  denotes  ordinary  absorption 
and  A^  extraordinary  (in  the  sense  of  the  magneto- ionic  components}: 


Table  I 

Night-Tine  Absorption  in  db  for  Several  Frequencies 


(mc/s) 

A,-  (  S  ?7Z) 
ca  ^  ju; 

Ac(>  375) 

Total 

AxU  375) 

Ay( >  375) 

Total 

SC 

C.G47 

.115 

.162 

40 

0.2 

c«4c 

.6( 

2C 

c.35 

1.9 

2.75 

0.91 

1.98 

2.89  ; 

14 

1.53 

3*9 

5-43 

1.8 

4.35 

6.15 

1C 

o  f- 
y 

7.4 

il-4 

4.2 

9.1 

13.3  | 

7 

7.2 

15.7 

22. 5- 

12.  C 

24.7 

3c.7  f 

Table  II 


Ordinary 

absorption 

■  »’  <1 
—  .  v*  — 

for  Several 

Frequencies 

f("'/  ■ 

s)  A:A<265}  A0(1C:-2CO) 

OUC> 

Total 

.fU(200-jlC) 

Ao(31C-700) 

Sub 

Total 

Total 

C.c.0 

C.16 

H..S- 

O.25 

C.17 

•42 

.64 

r 

1.01 

c,7- 

1 

J-«i  | 

1.C1 

Cc55 

1.60 

3*37 

20 

5,78 

3.02 

c.30 

4.34 

3.33 

7.6? 

14.47 

X4 

7.37 

6,c6 

13.43 

10.08 

7.7 

17*78 

31.21 

10 

14.3 

12.2 

26.5 

42.5 

15.7 

58.2 

84.7 

t«c- 


ii'/rn/  bt>“T 


Figure  3 
Dajr-Tia* 

Ordinary  asd  Kxteacrdinaztf 

Absorption,  t».  Frj^ueiicy 
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/igure  6 
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Table  III 

Dsy-Titce  Extraordinary  Absorption  in  db  for  Several  Frequencies 


f(gs-) 

*x(  >105) 

^(103-200) 

Ax(200-310) 

Ax(3io-?oo) 

Total 

80 

.26 

.19 

.25 

,20 

.90 

40 

1,03 

,80 

1.06 

.79 

3*70 

20 

4*08 

302 

4,72 

3 ‘53 

13  ‘70 

14 

8.19 

6.18 

11.63 

8.19 

34.19 

The  frequency  dependence  of  absorption  may  be  demonstrated  visually 
by  plots  of  the  data  presented  in  the  Tables.  Also,  plots  for 
differential  absorption  have  been  included  (See  Figures  3,  4*  5»  &)• 

The  determination  of  polynotninal  relations  is  relatively  easy , 
though  labourious,  once  given  a  set  of  points,  but  these  formulas  would 
obviously  be  merely  a  best  fit.  for  the  data  and  would  have  to  be  changed 
any  time  a  new  point  is  added  to  the  data,  so  for  the  curves  given  in 
the  log-log  graphs,  no  attempt  has  been  made  of  doing  this* 

Yet,  it  can  be  seen  that  the  plot  for  total  absorption  for  the 
ordinary  ray  at  night  is  a  fairly  straight  line.  So  that  ^  ^ • 
A  determination  of  this  leads  to  >7  =  2,0.  For  the  extraordinary  ray, 
the  first  three  points  form  also  approximately  a  straight  line  with  a 
slopfof  «=  1,4*  For  differential  absorption,  although  there  is  more 
spread  in  the  points,  a  straight  line  can  also  be  drawn  to  connect  the 
points.  The  slope  of  tnxs  line  is  v>  =-  4«3*  *or  day-time  conditions, 
no  such  determination  is  possib^  for  the  curve  joining  the  points  has 
a  slight  curvature.  In  this  case,  A  is  a  more  complicated  function  of 
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frequency  than  /-  7^  . 

Therefore,  for  night-time  conditions  as  assumed  here,  At~  ~~} 

J 

and  Aj>  ■?  ~ ,  in.  the  most  general  forms.  That  the  exponents  are 

f 

different  was  to  be  expected  since  differential  absorption  goes  to 
o,  or  at  least  approaches  undetectable  differences,  when  the  frequency 
increases,  Tnis  is  confirmed  by  the  fact  that  the  slope  determined 
for  differential  absorption  was  A  -  4»3*  which  means  that  it  decreases 
faster  than  either  ordinary  or  extraordinary  absorption. 

The  situation  is  somewhat  more  complicated  under  day-time  conditions. 
The  curve  for  ordinary  and  extraordinary  absorption  are  noc  straight 
lines  because  of  the  deviation  of  the  lower  points  from  a  straight  line# 
Tnis  is  due,  no  doubt,  to  the  presence  of  deviative  absorption  which 
for  longitudinal  propagation  does  not  follow  the  inverse  square  law 
for  the  frequency  and  would  be  expected  to  depart  unuer  general  conditions 
from  any  such  simple  behaviour.  In  fact,  the  upper  points  of  ordinary 
absorption  for  day-time  conditions  also  follow  an  inverse  square  law. 

The  differential  absorption  follows  a  straight  line  despite  some  scatter 
in  points  vita  gn  si  =  2.29  which  shows  a  slower  decrease  with  frequency 
tcan  that  occur ing  under  night-time  conditions. 

It  should  also  be  notea  that  the  ratio  of  the  absorption  at  night 
below  37 5  icm  to  that  over  375  fcsa  is  fairly  constant  for  all  frequencies 
for  both  the  ordinary  and  the  extraordinary  ray  (See  figure  7)«  This 
ratio  is  about  .44*  It  was  expected  that  it  would  be  constant,  given 


tas  form  for  absorption  4 


1 1 


It  is  also  seen  that  this  ratio  is 
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about  the  same  as  the  ratio  of  the  electron  content  below  375  to 
that  above  375  km,  an  observation  which  is  confirmed  by  Faraday 
rotation  determinations  of  the  electron  content  in  the  ionosphere. 

D  -  On  the  Height  Profiles  of  M  and  7 

The  day-time  attenuation  curves  show  that  there  are  three  main 
absorption  components,  one  at  about  3°0  km  and  the  other  two  at 
approximately  80  km  and  105  km.  The  upper  peak  corresponds  to  F -region 
absorption,  the  other  two  to  D  and  E  region  absorption.  These  components 
were  detected  by  i4itra  and  Shain  (Ref  35)  in  measurements  of  ionospheric 
absorption  at  18. 3  mc/s,  although  they  separate  these  components  into 
only  two:  F-regj.on  absorption  and  the  rest  of  the  ao sorption  which 
they  say  is  mainly  D-region  absorption.  The  results  presented  in 
Appendix  G  and  Table  II  suggest  that  E-region  absorption  is  also  important. 

For  night-time  absorption,  there  is  only  a  broad  maximum  of  absorp¬ 
tion  corresponding  to  the  F-layer. 

At  high  altitudes  the  difference  between  ^  and  ^  becomes  nearly 
constant.  It  was  shown  in  previous  chapters,  that  for  longitudinal 
propagation  and  high  frequencies,  this  was  to  be  expected.  The  fact 
that  it  also  occurs  for  a  case  of  general  propagation  suggests  that  it 
has  a  more  general  range  of  validity.  If  fcjj  ft  -  /cj  is  constant, 
then  fA  -  or  f^x  =  *  so  that  for  ranges  where  this  has  been 

determined  to  be  true  only  one  integration  need  be  performed  to  calculate 
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At  low  frequencies  (less  than  the  critical  frequency  of  the  layer) 
the  behaviour  of^  and  ^  is  very  complicated.  Absorption  for  these 
frequencies  could  also  be  determined  but  they  require  more  advanced 
methods  of  computation  based  on  contour  integrals,  (See  Sudden,  Ch,  20), 

The^  profiles  suggest  that  under  certain  conditions  waves  may 
be  trapped  at  a  certain  height.  For  example,  at  the  conditions  assumed 
in  this  paper  for  day-time,  a  wave  of  frequency  10  rac/s  and  ordinary 
polarization,  if  it  orginates  at  an  altitude  of  about  3°°  tan  on  an 
oblique  path,  will,  on  going  to  higher  altitudes,  find  a  higher  value 
of  ^  and  hence  will  be  refracted  down,  on  traveling  down,  however, 
it  will  again  encounter  higher  values  of  ^  and  hence  w'ill  be  refracted 
upward.  It  will  stay  thus  within  the  limits  imposed  by  the  ^  profile 
until  it  is  absorbed. 

Thus  by  these  brief  comments,  one  may  sse  that  the  height  profiles 
of  and  alone  merit  some  interest.  These  curves  have  been 
investigated  extensively  as  a  function  of  /V  and  ,  not  of  k  ,  in 
Buddea*s  and  Batcliffe's  book. 

It  may  be  added  that  the  profiles  presented  in  this  paper  were 
those  competed  by  using  the  Sen-Wyller  magneto-ionic  formula,  but  that 
comparison  computations  were  made  using  the  classical  Appleton-Hartree 
formula.  It  was  mentioned  previously  that  Shkarofilcy  had  showed  that 
the  classical  Appleton-Hartree  equation  could  be  died  in  the  high 
ionosphere  provided  one  used  the  appropriate  collision  frequency.  It  is 
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also  shown  in  the  Sen-Wyller  paper  that  in  the  limit  )7<.icO  the 
Ap>  oton-Hartree  could  be  retained  provided  one  used  a  collision  V 
frequency  equal  to  5/2  ^  (that  used  in  the  Sen-Wyller 

formula)  and  in  the  limit  the  appleton-Hartree  could  be 

retained  if  ^  was  set  equal  to  3/2  K.  ,  It  was  found  that 
\s*/  ~  for  tlle  kei£ht  range,  where  and  are  the 

absorption  coefficients  as  computed  by  the  Sen-Wyller  and  the  Appleton- 
Hartree  formulations  respectively.  If  the  Sen-Wyller  coefficients  are 
used  to  compute  the  total  absorption  then  a  certain  error  results  which, 
using  Shkarofsky's  remarks  and  the  results  mentioned  above,  may  be  easily 
seen  to  be  JEz /S'  (  *  where  /  <r,  .  c/A  is  the  absorption  as 

computed  from  the  Appleton-Hartree  formulation  in  the  region, 

which  may  be  assumed  to  be  from  $10  km  up.  That  is,  the  1  error 

would  be  Ji  where  is  the  £  region  absorption 

as  computed  from  the  Sen-Wyller  formulation.  This  error  is  an  over- 
estimation  of  the  absorption  since  in  the  £  region  absorption  coefficients 
which  are  2.5  higher  than  the  more  appropriate  ones,  then  /*j  fAf£)  -  /ej  4 
represents  the  shift  in  log  coordinates  where  A  is  the  correct 
absorption,  or  P-  t °J  Cp~ £ )  where  Z5  is  the  plotted  absorption. 

3ut  y?"  /0J  (/-  JLJ  where  «  f  .  Hence 

/«*  P  -  J'j(r-E)  -  "  tf)  ~  -f 

where  terms  of  order  two  or  higher  have  been  neglected.  This  ratio  is 


Hence 


('?)  -  f 


not  only  small,  .13  for  extraordinary  day-time  absorption  and  .12  for 
tne  ordinary  on  the  average,  but,  perhaps  surprisingly  independent  of 
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frequency,  hence,  the  expressions  determined  for  frequency  dependence 
are  not  affected  because  the  shift  is  constant  for  each  ordinate.  It 
there*’ ore,  appears  to  justify  the  use  of  Sen-Wyller  only,  although,  it 
is  realized  that  in  some  other  future  study  both  should  be  used. 

The  fact  that  =  2*5  the  present  computations  is  by  no 

means  a  general  relationship  but  the  conditions  under  which  it  is 
obtained  should  be  investigated.  Time  precluded  its  investigation  in 
the  present  paper.  Because  in  most  of  the  ionosphere  y<*0  for  the 
present  computations,  then  if  =  S ^  as  mentioned  previously, 

the  Appleton-Hartree  equation  can  be  retained;  in  general  cases,  it  is 
not  readily  possible  to  say  what  this  implies  because  of  the  complicated 
algebraic  formula  but  in  a  simple  and  useful  case  one  can  derive  certain 
requirements,  i.e.  the  longitudinal  case:-  in  that  case  and  for  , 

f^Ati  written  in  terms  of  the  frequencies  is: 

•>t<) iCu+s)*-*  *AH] 


while 


then  if  y  -  Jf  uL 


| 


I 
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c  n.\4S  j 

30  that  i.y^  (  -y-1  J  "u  ^  within  the  limits  of  error  obtained 
by  Sen  and  Wyller  (Kef  49 )«  If  were  x?^  as  was  assumed  for 
the  Apple ton-Hartree  computations  that  were  made,  then 


and  %  \a«  -  ^ 


so  that  - 

propagation  that 


then 

7^7 

would  imply  for  longitudinal 
/ 


E  -  On  /C  or  Absorption  Per  Electron 

Some  aspects  of  K  may  be  explained  qualitatively  if  use  is  made 
of  several  results  developed  in  Chapter  3,  Equation  (3-24)  gave  £  as  a 
function  of  frequency  and  it  was  shown  that  K  has  a  maximum  at  » 

and  it  was  therefore  explained  that  this  result  implied  that  the  maximum 
for  K  vaa  higher  for  the  extraordinary  ray  than  for  the  ordinary.  The 
graphs  in  Appendix  H  show,  the  same  type  of  variation.  Tnat  is,  their  maximum 
is  frequency  dependent,  and  the  higher  tho  frequency  the  lower  the  maximum i 
also  Kx  is  generally  higher  than  K0  although  for  low  frequencies  the 
behaviour  of  both  of  these  becomes  erratic  and  they  cross  over. 

The  main  question  asked  about  was  whether  it-  was  proportional  to 
/■/or  not,  i.e.  whether  ^  -  KN where  *C  depends  only  on  the  frequency  and 
the  collision  frequency  profile.  The  computed  results  show  that  for  the 
Sen-Wyller  approach  (and  hence  for  the  Appleton-Hartree  formulation 
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according  to  the  comments  made  before)  and  under  the  conditions  assumed 
in  this  specific  problem  (dip  angle  of  73°  <*Qd  gyrofrequency  of  1.57  mc/s 
corresponding  to  Sagamore  Hill  Observatory,  Mass*  and  the  assumed  profiles) 
the  linearity  'ition  for  the  ordinary  ray  holds  vithin  certain  limits. 

It  is  good  eiiv.u-.gh  for  frequencies  above  2  mc/s  at  low  heights  where  the 
collision  frequency  is  ^igh  and  the  electron  density  is  .low.  As  /V 
increases  the  linearity  condition  becomes  less  valid.  It  may  be  said 
that  the  linearity  condition  may  be  taken  to  be-  valid  up  to  about 
for  frequencies  down  to  2  mc/s  and  up  to  10^*  for  frequencies  down  to 
10  mc/s.  ?or  high  frequencies  such  a.s  10,  4-»  ard  80  mc/s  it  may  be 
taken  to  be  valid  at  densities  of  10®  /cc  which  are  densities  which 
occur  in  Polar  Cap  Absorption  events.  Outside  these  limits  deviations 
from  linearity  may  be  considerable* 

As  a  side  result,  it  may  be  noted  that  the  height  at  which 
differential  absorption  appears  is  dependent  on  the  collision  frequency 
and  the  frequency  of  the  wave,  so  that  if  a  rocket  probe  with  a  transmitter 
at  a  certain  frequency  wave  sent  into  the  ionosphere,  then  that  hei^it 
at  which  the  powers  received  in  the  two  wave  modes  begin  to  differ 
appreciably  is  characteristic  of  a  certain  collision  frequency.  If  the 
transmitter  were  a  multifrequency  one,  some  mapping  of  the  collision 
frequency  profile  at  low  heights  could  be  done.  The  idea  is  worth 
considering* 


A  Method  for  Measuring  Cosmic  iioise  Levels  and  Total  Absorption 

It  tuxs  been  mentioned  before  that  the  4^  -  foj  /p0  was  largely 


104 


SP/PH/65-1 


especially,  in  the  high  ionosphere  as  was  shown  by  the  plots  in  Appendix 


G.  This  implies  that. 


since 


,  then 


Where  the  constant  L,  depends  only  on  frequency.  Therefore,  this 
constant  once  theoretically  known  permits  the  assertion  that  P~ 
where  ^(J )  i3  a  percentage.  Now,  since  D  is  measurable  by  riometers, 
then  dp  can  be  direct ly  calculated  and  hence  the  cosmic  noise  level.  By 
way  of  illustration,  if  ''/'(J)-  »1  and  it  is  determined  that  $■=.  /  pB  , 
then  is  10  b B  and  from  A0  -  P~  ,  where  P  is  the  cosmic  noise 
level  and  ^  the  power  receiver  in  the  ordinary  component  above  a  certain 
reference  level,  then  P  =  /o  t  ?0  .  This  method  is  close  to  that  of 

Little,  etal  (Ref  ).  Although  they  determined  by  extrapolation  of 

a  plot  of  %  vj.  .  For  more  details,  the  original  paper  should  be 

consulted. 

In  the  Appleton-Hartree  formulation  and  for  longitudinal  propagation e 

/  iJf,  )  2. 

for  aou-deviative  absorption,  ^  r  (  ,  so  that  it  is  very 

simple  to  compute  }(f)  for  any  frequency.  It  is  of  interest  to  see 
whether  propagation  in  any  general  direction  would  alter  this  situation 
very  much.  In  the  case  presented  here.  Appendix  L  and  A  tend  to  indicate 
that  it  does  not,  for  graph  3  in  Appendix  K,  indicate  that  for  a  given  N/  , 
tae  differential  absorption  once  it  appears  is  largely  insensitive  to  the 
collisi  requeccy  since  it  remains  constant  for  the  heights  examined. 
Appendix  G,  as  mentioned  before,  shcvs  Px  -  for  the  high 
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ionosphere  is  constant.  Lastly,  so  does  the  linearity  of  ^  at  low 
heights.  This  fact  makes  the  selection  of  the  profiles  fiXplained  in 
Chapter  IV  less  arbitrary  than  would  seem. 
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VI.  CONCLUSIONS  AND  KSC«-M£NDaTIONS 
The  main  point  that  this  paper  tries  to  make  is  that  differential 
absorption  may  serve  to  calculate  total  absorption  and  cosmic  noise 
levels.  For  example,  differential  absorption  is  directly  measurable 


and  if  theoretically  one  estimates  that  for  a  certain  frequency  the 
differential  absorption  repi’esents  10%  of  the  extraordinary  absorption 
then  for,  say,  /  differential  absorption,  the  extraordinary  absorp¬ 
tion  for  such  a  frequency  is  it  t&  ,  then  as  mentioned  previously,  the 
cosmic  noise  level  is  the  aum  of  /oiP>  and  the  power  in  AS  above  a 
certain  reference  level  received  in  the  extraordinary  mode. 

Also,  the  frequency  dependence  of  absorption  A ,  and  differential 

C 

absorption  was  determined  to  bas  yf  y,  «  j-i-3  ,  respectively  at  night¬ 

time.  during  day-time  conditions,  no  formula  appears  to  fit  the  data. 


but  the  upper  points  of  the  ordinary  absorption  curve  appears  to  obey 

C*-  v  f 

the  law,  »  differential  absorption  was  proportional  to  ya-*  • 

The  approximation  j-r/r At  used  often  in  ionospheric  physics  is  shewed 
to  be  valid  under  general  conditions  at  low  heights  and  high  frequencies 
and  for  low  //  «  Outside  those  limits  considerable  departures  can  occur# 


rsecasmendations : 

It  is  recommended  that  a  systematic  study  be  undertaken  of  the 

x  «*'  M  o  ,  _ 

variations  undergone  by  _ as  the  angle  of  propagation  varies  for 

different  typical  profiles,  with  the  intention  of  telling  how  much  independence 
from  the  profiles  actually  exist  and  how  they  affect  the  percentage  differential 
absorption.  Furthermore,  the  absorption  calculations  should  be  combined 
with  a  ray  path  analysis  to  take  into  account  deviative  losses  with  more 
accuracy. 
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£a  experiment  should  be  performed  after  better  calculations  are 

Ax  -A, 

done  of  — -- —  to  determine  total  absorption  from  percent 

r  a 

differential  absorption  and  also  cosmic  noise  levels  outside  the 


to  determine  total  absorption  from  percent 


ionosphere.  Lastly,  rocket  measurements  as  low  frequencies  could  be 


used  to  determine  collision  frequency  profiles. 
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APPEcZDEC  A 


The  Wave  Magnetic  Field 

wave  traveling  in  the  Z  direction  and  linearly  polarized  with 
its  electric  vector  in  the  x-direction»  will  produce  a  force  on  a»'. 
electron  in  the  x-direction 


m  ux 


IK  s  Cosuft 


(A.1) 


r  sin  <0~t 


(A.2) 


frco  equation  (1.15),  and  taking  j(  •  where  M  is  the 

refractive  index  (in  general,  complex)  and  n  is  a  unit  vector  in  the 
direction  of  propagation,  2.  ,  ~  *  ,  c  is  the  velocity  of  light, 

the  following  is  obtained:: 


ME*  -  c\ ff 


M  -  HftL 

7  /^r 

Therefore,  the  magnetic  force  experienced  by  the  electron  in  the  direction  y 

<r  *  tM.  (~j) c,i  u' 

~  ^  s in  n't c»s  tit 

hnuPC 
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This  force  is  a  maximum  when  S/» 


F  .  *f.*M 

rM  -  - - 

iM  COG 


The  ratio  of  the  magnetic  force  to  the  maximum  electric  force,  is: 

f”/f<  TJGZTc/'*' 

,  t£.M 

i  ht tOC 

"The  value  of  £r,  at  100  Km  from  a  transmitter  radiating  10^W  is  of  the 
order  0.08  lM%  and  the  magnetic  field  of  the  wave  is  about  2.7  x  10**^ 
oersted.  The  magnetic  and  electric  fields  encountered  in  the  ionosphere 
from  man-made  radio  transmitters  will  rarely  reach  these  values,  and  the 
frequency  is  always  greate'  10  kc/s.  The  refractive  index  is  of  the 
order  of  unity.  Hence  the  ratio  will  not  exceed  3*7  x  10*"^",  As  can 
be  seen  this  is  a  small  ratio  in  the  "worst"  of  circumstances.  This 
justifies  it  being  neglected  in  the  equations  of  notion. 
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APPENDIX  B 


I-li-i  >i "  >  1  )T-\5  W  >1 


B!i«Sv5iniE*33U 


The  dispersion  .relation  (ill)  in  the  Apple  ton -Har  tree  formulation 


S/HXS  Sih  *±  +  ZCCJXfi 

J--J.  J--JL-  JL  -JL  '  ° 
*  m 4  **  **  /** 


(111) 


where 


a  x _ _ 

£,  r  /  -  r  /-  ~r~7~~7 

p+ji  t+  r-tZ 


(B.l) 


**  5  / " 


r  /  - 


/-K-/Z 


(B.2) 


*  . _ JL 

*3  '  l~  j9 


Equation  (111)  may  be  recast  in  the  form 


(B.3.- 


(.«,»<*  (**''**•  .  » **  fgl*  _ 

*  ^  **  —  r"  — '  — 

Mv- M*-  *y 


The  summation  of  these  leads  to 


((, + (>)J+  2 (,  fx  f ,  -* 


(B.4) 


(B.5) 
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A  series  of  redefinitions  is  useful  at  this  point*; 


=  *1*2 


B  =  2g, 


=  «1(*3-i2)+  t2(a3-«1) 


D  =  BA 


E  =  2T 


K2*3 


(B.6) 

(B.t) 

(8,8) 

(B.9) 

(WO 


so  that  equation  (B,5)  then  becomes :  V 

+  MtCs^1»C-  tj  1-  £. j*  (B>U 

•  ’ 

* 

•n 

If  use  is  made  of  the  formula  for  determining  thfc  ;#olutions  for  & 
q.uadtatic  equation,  then  ’  : ; 

4  VtoS*C)+  C Si»X9-  >  +  S] 

*  "  *  2  (^siw1*  4- 8. c *>$*•)“ D  JT  J(c$inx9 ->)*- 4 (A 8 <«V)£ 

» 

The  terms  inside  the  radical  may  be  sean  to  be*  •  . 

.  ’’  -£v 

(s.ij) 


But  by  the  definitions’  of  A,  Bt  C,  D, 


and  2CD+4AE-4B*  may  be  seen  to  be  4  (€*)*■  (t,-  **$!--  ^  '  tbis  last 

.ife?*:..  ■  -  -  , 

expression  is  set  s%uel  to  iP^,  90  that  the  terms  the  radios! • 


become 


V** 
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The  use  of  E^,  Ig  and  is  no  longer  useful,  so  from  here  on  the 
ionospheric  parameters  X,  Y,  Z  are  used.  With  this  in  mind. 


7  /  /-«•*/  lo-r-iz)(i*Y-iz)l  rj 


Also 


[lJLIL _ - — —  / 


(B.17) 


So  that  the  terms  inside  the  radical  divided  by  4  can  be  expressed 


n-x-izl'rjJZ _ tl  Yr*  .Wl 

[TTIz  J  U('-x-ar  “> 


0*1$) 


The  term  in  the  numerator  of  the  fraction  may  be  written  ass 


rxaJ 


(B.19) 


and  the  other  term  in  the  denominator  is 

r  xY1  ,  zxo-iz-t) 

Ifi-Tajfcw wr* *  ~Y-aUY-^) 


(B.20) 


So  that,  finally, 

lO-iZ'XJlL 


1  -  u-Lz~YXi-a*vX>-iZ} 


(£.21) 


U7 
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It  follows  immediately  from  (B.  21)  that 


V _ 

lO-X-cZ) 


i/2 

I- 


This  is  the  Appleton-Hartree  equation# 
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APPENDIX  C 

Determination  of  P*..  K\  %,  ^  Sx 

Boltzmann's  equation  becomes,  when  reduced  by  the  simplifications 
made  in  Chapter  II:: 

4  [>,  vj,  --  fff  (f.v;-  a  jjUJcJa , 


Furthermore  in  the  Chapman-Enskog  method  of  solution 
fi  =  ^4.  Pa  *Va  (**  CQZltik+piStnuSk  )-f  SUUtjt  J 

+*  L’H0x(H.xPa)].Vfc(lVcoswt4  ) 


(2.7) 


In  principle  one  substitutes  equation  (2.7)  in  Boltzmann's  equation 

and  equates  the  coefficients  of  p  *  fV .  V*.  »  C  H.»  A  P»  )  •  and 

)  *-  — *  —  — 

{Jtf.  *(&  xPaJ]  -  Vx  on  both  sides  of  the  equation  to  each  other.  It 

is  easier  to  do  it  term  by  term  instead  of  just  outright  substitution. 
It  should  be  remembered  that  for  any  function  f(r)f 

*  Ur)  =  if  V  r 

r  \a*  wr  • 


?r 


(c.D 


Therefore 


•  £  (f) 


(C.2) 


as  a  help  in  differentiating  dots  products  of  V*  with  some  other  vector* 
7  f,  can  be  written  also  as  rrr  which  merely  means  that  jiT 

V  *-  #TJ| 

is  the  vector  operator  whose  i'll*  component  is  ,  so  that  when  it 
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operates  on  a  dot  product  where  the  j'-s  imply  summation, 

then  (A;  U1-)  -  A;  «  With  these  considerations  in  mind  then 

*  i 

}  j  :  M  >Jl‘\  l\.^(£coiUt  +  ifcsihvtyLfcayuUp'SiMVl 

‘  tf  {*4  ~  ~  iir>  ! 

+  +Qsi «"%.<!»  v,'  f 

4-_^|^*coJwt  +  —  si«“f  j  |H.  #  (M fj)J '  *i}+  ( t  t0  3 


+- 


ihtoUQh 


(0.3) 


Therefore 


£>C'Sutv*j,  =  |i  f  l  K(f£tosJ+i^sj»ut) 

+  (^osarf  +  Q*i»wt)(S,*[‘)-*x  +(£s'iU>i  +  £s'”“>V 

Xf(^JX(«.<£i).VA|+  f^CcstdC^toseot  i p^Sincot)  - 


r2  coswf  +  S2si>uot)n*fl**ih3-iJ’  “Siat 


(0.4) 


and  also 

^  -  5;  /  =  >A  SfHtitt**  x(j.)  ■  /*  *. 

>>7^  / 

^  (f*  c<?W s‘'»e*tJlh.x(&*£x)]'V9.  - 

tC (*zf$ujt+  K9'»«ft)/{(t.x£)-  V^J 

(c.5) 
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and  finally 

—  -  i*  '  K  +  LUf3x  COStof)  +.  (t!„  Xf\)  'j£  (-M  fi  Stftiot  + 

io»L  tostoij + Of*  *(M>  [(-  &  /i  s/> 7  ^  z9  + 


UJ$X  i.OS£vt ]  Co.  6) 

So  that  to  pick  the  coefficients  of  l\  *  / 1  *  J4  , 

and  i id»  *(&»X {*)]  'V*.  one  has  to  look  in  (C.4)»  (C.5),  and 

(C.6).  In  this  fashion  the  equations  obtained  are: 

Pz*\/J cosrtt^i  \  U)  LI  (-<4  sinui +8^  cos  >*+)  l  = 
f[(  [{/(*!  COStftf ^ si* ut )(£  .)&'  )- 

f  (c^cosu^  f  jixsi»0t  i(q ,  va)J  j  y  We  dv, 


-  j-,  (i^ssui  4-  ^smut)[(>(.%i\)-v,J  I  (0>8J 

&•  <  ( H  •  x  Pz)]  •  Vi  j  ^  c  %».  cos  coi  +  +to( -t^irudt  +iac©*u$t)| 

-  ///  {I'cr^osiot^.  St  s;»ui)[(!t.x(!!*  *£))•  vj]  — 

t  (/l  c *s  u>t + £  si  tut) [Of,  x(& x$J)-V>Jjl,Jl<f€  dv,  (c.9) 


and 

l?*Ci>Slct[(<**COSiOt+J5zsinClrf'}  +  cost ot  4  IfUitHtoi ) 

-  HtstoVjOi  COSKit +i%si»ui]  /  ^fAcesa)1>  ¥ 

'  »cra 

£rSi„at)l]=  iffa/C-  /.C'jjiSM  (c.xo) 
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The  formidable  integrals  on  the  right  have  been  evaluated  by 
Chapman  and  Cowling  (Refi*  )  using  several  approximations.  First 
of  all,  )4  " V2/ ^  j  (the  relative  velocity).  This  is  equivalent  to 
the  assumption  of  infinite  mass  for  the  heavy  neutral  particles 
(anologous  to  Rutherford  scattering).  Secondly,  f ,  is  assumed  to  be 
Maxwellian  as  was  pointed  out  in  Ch.  II,  Section  A.  Thxs  considerably 
simplifies  the  integrals;  for  example,  the  one  in  (C.7)  reduces  to 


(H  COS  Oft + j3x  sintdix 71  'J4  jfff '(£  '■  {,}  f  b  (tUe  A  (C .  11) 

The  last  integral  is  equal  to  according  to 

Chapman  and  Cowling.  (C.S)  and  (C.9)  can  be  treated  the  same  way. 
In  equation  (C.10)  a  different  procedure  is  followed::  each  side  is 
multiplied  by  Jv*.  c  4%  2 ;  the  resulting  integral  on 

the  right  has  also  been  evaluated  by  Chapman  and  Cowling  and  is 


(C.12) 


All  these  equations  lead  to  a  system  of 
The  system  ist 


linear 


equations  for«y^  ^ 


(C.X3) 


(c.14) 


CC.X 5) 
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.nir’T'Mrvrr  n 


The  Scn->-/yller  formula 

In  Chapter  II  the  dispersion  relation  for  the  index  of  refraction 
of  the  medium  v;as  given  as: 

[-(*!.  +  %  ^rr  CoS4 


-  („  C0S<p 


Ml-(*r+<3?)  errS'^t 


Sitij  ces$  —  5  *  ^  4  M- (£j  +  %Sl>,i^) 


=  0 


(2.57) 


Evaluation  of  the  determinant  leads  to: 

C°s<f{[>IIcw$XM  “  (ti+tfU  4-  %^{TT  Si»X$Z°s<}  [  + 


(m  ^cos4{4s^4 =  0 


Algebraic  Manipulation  of  the  previous  equation  leads  tot 

„*(-2-£us^)+  M1^  +  ac‘,jl+)  '  |  fB  +  f]  a  0 

Where  .a,  B,  C,  D ,  E  are  as  defined  in  Chapter  II,  If  the  quadratic 
equation  for  is  solved  by  the  usual  formula,  then 

_/yi4lscori4)t  jAt+2*'&coslb+tf’cos*4 -  Q 4 

-(>+•£  cos*-$ ) 

_/Pte  ) 

If  the  fact-  that  -i4B  -i&sK^T'  /  =  0  is  used  along  with  the 

equality  CX:  then  the  quantity  under  the  radical 

reduces  to:; 

S*Co5^  -  C 
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Therefore  _ _____ 

,  A+Bcosx4>  +  Jb1coi46  - 

M  r  — >— - — -  -  f _ 

D+-  £  cos1 4 

And  since  the  angle  of  rotation  <f)  *  <f-  ^  ,  where  (f  is  the  angle 
between  the  propagation  vector  and  the  magnetic  field,  then 

i  A -t zcp  +  J IS** '»  V  -  Czcos*a~ 

^  T>+  E  s/nx<fi 

This  is  the  Sen-Wyller  magneto-ionic  foraula. 
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APPENDIX  E 

A.  Dipole  Field  Deviation 

It  is  desired  to  prove  that  for  a  dipole  field  the  angle  between 
the  radial  direction  and  the  lines  of  force  is  constant. 

If  a  spherical  coordinate  system  is  used  with  origin  at  the  center 
of  the  earth  then  a  line  of  force  for  a  dipole  field  is  given  by 


t-  l  s/>t  t-6 


(E.l) 


Where  jr  is  the  distance  from  the  center  of  the  earth,  Q  is  the 
colatitude  and  /  is  a. parameter  which  characterizes  each  line  of  force. 
A  function  ^  is  defined  such  that 


^  r-  L 


(E.2) 


Since  the  gradient  of  a  function  in  a.  spherical  coordinate  system 
ia  given  by 


3$  t  ilt  1  iJL  I 

V?  -  ( Jr~t  r"  J# ,  rsi»6  ' 


(E  .3) 


Then 


Vf  *  0,  'ZCoto>  0 )s  (l,  *  }  (E.4) 

or  using  (E.l) 


7^  =■  (t  -  2  cot 9^  o  ) 

The  vector  f  f  is  normal  to  ^  =  constant,  in  particular,  it  is  normal 
to  ^  =  0,  that  is,  normal  to  the  lines  of  force. 

For  spherical  surfaces,  i.e. 


V<?  -  (\}0}  0) 


(E.5) 
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Hence  vf-  VF  =1  =  a  constant.  Any  vector  A  parallel  to  the 
lines  of  force  at  a  point  satisfies  the  condition  A  ■  vt-  o 
In  particular,  A  can  be  the  magnetic  field,  since 


„  ,  /-iMctsi  -Ms  in e  I  /  . 

a-vf  --  ( — — ,  —  oJ-0.  -kute‘  °> 

j  '  j 

-  --2 MtesB  ^XicosG 

r  ~  r*  r* 

H-V'+cO 

Therefore,  fl  x  (?rx  J  —(£■  Vjt) vr-  (#- ?r)vt 
or  by  (E.7) 

AxfrrxVt)  -  -(tt.vr)yt 

but  Hr  and  VrXVft  =  (  6  Cf  -  i  cot-0  )« 


Therefore 


(Hrj  H0J  o)K(Q ;  <7y  J  cite )  s  Hr 


(E.6) 


(E.7) 


(E.8) 


(E.9) 


Ucoi*#9j  2  Hr  tote )  o)  -  tffVt 


(E. 10) 


Two  equations  result 


2toteff9  =  -fi  Mr 


(E.n) 


Hr  -  Hr 


(E.12) 


The  second  one  is  trivial,  but  from  the  first  one,  one  can  see  that 
since  Hr  -  where  £  is  the  angle  that  H  makes  with  tl 


radial  direction  then 


t  - 

>  ‘  H 


JL  n 

2  co  tin  o 


(E.13) 


(3.12j) 


jrSaJSgjjgsg 
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COS 


J 


- 1  CoJ$ 
(f-t-iC6SZ0)  V* 


from  which  it  is  evident  that  £  is  independent  of  which  was  to 

be  proved. 
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APPENDIX  g 

Individual  Height  Profiles  of  M,  *L 

Because  they  are  so  numerous,  it  was  decided  to  place  the  results 
in  an  Appendix*  The  figures  are  almost  self-explanatory  once  it  is 
’acwn  what  the  symbol?  stand  for,  A  stands  for  electron  density,  B 
for  collision  frequency,  0,  A  for  and  ^  respectively  while  + 
and  -  denote  ordinary  and  extraordinary  ray*  respectively.  The 
frequencies  and  the  conditions  (day-time  or  night-  )  for  which  the 
profile  is  computed  are  set  down  at  the  top  of  the  graph.  The  ordinate 
in  all  of  these  is  the  height.  The  graphs  are  arranged  in  blocks,  one 
for  each  abscissa,  and  written  each  in  a  progression  of  frequencies  from 
the  lower  .to  the  higher  values. 
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APKuNDIX  Qr 

Grouped  Height  Profiles  of  U.  Z 

Part  of  the  graphs  presented  in  the  previous  appendix  are  presented  but 
grouped  together  according  to  their  being  JQpt  3$o.  Only  those  which 
were  used  in  the  determination  of  frequency  dependence  are  plotted  ^re. 
These  aro  clearly  self-explanatory*  The  grouping  makes  fcbeir  frequency 
variation  easier  to  see. 
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iiP  j  j-l  iDIX  H 

This  appendix  contains  various  graphs  examining  the  relation  between 
N  and  i.e.  under  what  conditions  linearity  holds.  The  frequencies  used 
are  stated  in  each  of  the  graphs.  The  first  two  present  20(log^Qe)lC  for 
both  the  extraordinary  and  the  ordinary  ray  against  hei#it.  They  are  pre¬ 
sented  for  both  almost  longitudinal  conditions  (THETA  =  .02  radians)  and 
for  almost-  transverse  propagation  (THETA  =  1.56  radians).  The  third  graph 
shows  for  both  l(x,  l^o,  their  variation  with  height  when  6  =  1.27  and  N  =  1. 
The  angle  6  =  1,2?  was  the  one  used  in  the  general  calculations  presented 
before .  The  last  one  shows  a  comparison  (for  ©  =  1*27)  between  the  attenua- 
tions  given  when  N  -  1  and  when  N  =  10  /cc. 
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